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FOREWORD 


The work described in this report was performed by Information & Control Systems, 
Incorporated (ICS) under Contract Number NAS1-19091 for the National Aeronautics and 
Space Administration (NASA), Langley Research Center (LaRC). The work was sponsored 
under the Small Business Innovation Research (SBIR) program by the Aircraft Guidance 
& Controls Branch. 

Mr. Aaron Ostroff was the NASA Technical Representative for this contract. Dr. 
Nesim Halyo was the ICS project manager. 

All the design, analysis and simulation results shown in this report were obtained 
using our Computer-Aided Control System Design (CACSD) software package ACET™. 

We would like to thank the National Aeronautics and Space Administration for its 
support of this important work. 
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SUMMARY 


This work was performed under the NASA SBIR program. ICS plans to market the 
software developed as a new module in its commercial CACSD software package: ACET. 
This module will contain algorithms for designing SOFFT control laws, for analyzing com- 
mand tracking performance with Z-plots and for simulating SOFFT control systems using 
the incremental implementation. ICS also plans to market the SOFFT methodology to de- 
sign flight control systems for high-performance aircraft and in other control applications. 

A new control design methodology is proposed: Stochastic Optimal Feedforward and 
Feedback Technology (SOFFT). Traditional design techniques optimize a single cost func- 
tion (which expresses the design objectives) to obtain both the feedforward^ and feedback 
control laws. Since the feedforward response must be fast while the feedback response 
must be relatively slow (to attenuate noise), combining these objectives into a single cost 
function produces conflicting demands; so that neither of the objectives is fully achieved. 
In the SOFFT approach, two cost functions are defined. The feedforward control law is 
designed to optimize one cost function, the feedback optimizes the other. By separating 
the design objectives and decoupling the feedforward and feedback design processes, both 
objectives can be achieved fully. 

The main feedforward design objective is to produce a desirable response in tracking 
input commands when no random noises and disturbances are present. In particular, a 
fast and smooth tracking response during the transient phase while performing difficult 
maneuvers is the goal of the feedforward law. The main feedback design objectives are 
to suppress sensor noise, accommodate plant disturbances and provide stability in the 
presence of plant modeling uncertainties. 

A new measure of command tracking performance, Z-pIots, was developed. By ana- 
lyzing these plots at off-nominal operating points, the sensitivity or robustness of the total 
system in tracking commands can be predicted. 2-plots provide an important tool for 
designing robust control systems. 

The Variable-Gain SOFFT methodology was used to design a flight control system 
for the F/A-18 aircraft. The Variable-Gain SOFFT controller can be used to extend the 
operating regime of the aircraft and to provide greater performance (flying/handling qual- 
ities) throughout the extended regime. The angle-of-attack (a) command system designed 
smoothly and quickly brings the aircraft to within 3 degrees of its physical (actuator) limit 
with an easy maneuver. A detailed nonlinear simulation of the aircraft and control system 
displays excellent command tracking performance. 

An important by-product of the SOFFT approach is that by removing the conflicting 
demands on the control system, both feedforward and feedback control laws can be designed 
with less effort and in considerably less time. Our main conclusion is that the concept of 
separating the feedforward and feedback objectives and decoupling the two design processes 
works well. 
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I. INTRODUCTION 


While the feedback control problem has been studied extensively (e.g. [l]-[l2]), the 
feedforward control problem has received less attention [13], [14], [8]. Stochastic opti- 
mization methods are often used for feedback control law design (e.g., Linear-Quadratic- 
Gaussian (LQG), Stochastic Output Feedback, etc.). However, stochastic optimization of 
feedforward control systems has not been exploited to its full potential. 

In this report, we propose a new control design methodology, namely the Stochastic 
Optimal Feedforward and Feedback Technique, or simply SOFFT. The work was per- 
formed within the context of the SBIR program, and describes the research and devel- 
opment activities undertaken. SOFFT is an Integrated Feedforward and Feedback (IFF) 
design approach. Unlike currently used techniques such as model following, in SOFFT the 
feedforward and feedback control laws are designed independently of one another. The 
feedforward and feedback control laws are then integrated using the SOFFT structure. 
In explicit model following, the feedforward and feedback control laws are obtained from 
the optimization of a single criterion in which the performances of the feedforward and 
feedback control laws are jointly evaluated. As a result, the optimized criterion results in 
a compromise between the performance of the feedforward and feedback control laws. 

In other cases, the designer tries to meet all the control objectives by appropriate 
design of the feedback control law. Generally, this results in the error feedback structure 
shown in Figure 1. Trying to achieve all the control objectives using feedback 
alone places conflicting demands on the control law, making it difficult and 
sometimes impossible to achieve all of the objectives. Some undesirable and un- 
necessary compromises may have to be made with time-consuming trial-and-error designs. 

Thus, both in the error feedback approach and in model following, the design is the 
result of compromising feedforward and feedback objectives. However, is such a compro- 
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Figure 1. Error feedback controller 

mise really necessary? The SOFFT approach decouples the feedforward and feedback 
control design process by separating the feedforward and feedback control objectives. In 
the SOFFT approach, a compromise between feedforward and feedback is not necessary. 
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Figure 2. Explicit model following control law structure 


In the explicit model following approach, the resulting system structure has a feed- 
forward control law as shown in Figure 2. It is desired that the system respond like the 
command model which is now shown explicitly. On the other hand, the error feedback form 
may implicitly contain a desired command model, but the system response is determined 
by the closed-loop system time constants and natural frequencies. 

Consider the model following problem 
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Xk+1 = 4>x Xk + r* U x k + Wx* 


( 1 ) 


Zk+1 = 4>zZk + T z U z k (2) 

k = H x Hz %k (3) 

J = J™, 2 (jyT 1) g E e * g e * + “«* R “•* « 

where E denotes the statistical expectation operator. 

An explicit model following problem is defined by the plant in (l) , the command model 
in (2), the error defined by (3) and the cost or objective function given by (4). Assuming 
that {w xk ,k > 0} and { u z k,k > 0} are independent gaussian white noise sequences results 
in an optimization problem. The optimal control law has the structure shown in Figure 2 
and is given by ([ll] pp. 548-549) 

= ~K X x k + K 2 z k (5) 

Note that the feedback and feedforward gain matrices K\ and K 2 axe not indepen- 
dently designed, but are highly interdependent. In fact, both depend on the weighting 
matrices Q and R. Note how conflicting demands are placed by Q and R. Increasing Q 
means that tracking the input commands is more important relative to the control effort 
expended to achieve the tracking performance. As long as the control penalty term R is 
present, there will be tracking error even for constant commands in steady state. 

Now suppose that a difficult maneuver requiring high control effort is commanded. As 
long as actuator limits axe not exceeded, we see no reason for the feedforward controller to 
sacrifice agility or tracking performance in order to minimize the control effort. We think 


3 



that the feedforward control should expend the necessary control effort and perform the 
maneuver with as much accuracy as possible. However, in the explicit model following 
approach, such a feedforward control would also result in a feedback control law with very 
high gains and very little filtering or disturbance robustness. It is difficult to make sure 
that both a robust feedback control and a fast-response feedforward control are obtained 
by minimizing a single cost function weighting the tracking error and control effort. 

While a particular control design problem may contain additional objectives, most 
designs try to achieve the following overall control objectives: 

1. Quick response to input (e.g., pilot) commands during fast maneuvers 

2. Desired overshoot and damping characteristics 

3. Attenuation of (low response to) sensor noises and high frequency disturbances 
(e.g., turbulence) 

4. Non-oscillatory response to large disturbances (e.g., high gusts, shear winds, etc.) 

5. Maintain the above tracking and disturbance characteristics despite uncertainties, 
variations and nonlinearities in the plant, actuators and sensors. 

The error feedback controller structure is unable to achieve all of the above 
objectives. In particular, objectives 1 and 3 (and often 5) place conflicting 
demands on the feedback controller. For example, suppose that the input command 
is a step or a pulse. Initially, the error e which drives the controller ( F ), will suddenly 
increase. In this scenario, the controller must let this sudden (high-frequency) command 
pass through so that a quick response to this command can be achieved (objective l). 

Now consider the scenario in which the command remains constant but the measure- 
ment noise has a random spike. The error, e, will again suddenly increase. However, 
in this case, the controller must not pass this signal through to the response; it must in 
fact attenuate it (objective 3). Thus, conflicting demands are being placed on the con- 
troller, F. Clearly, the error feedback controller cannot meet both objectives. Therefore, 
a compromise must be reached in which neither objective is completely achieved. 
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NOISE 

Figure 3. SOFFT integrated feedforward and feedback (IFF) control structure 


In the SOFFT approach, the control objectives are separated into two 
groups. Objectives related to the system response to tracking commands 
are met by the feedforward controller. Objectives related to disturbance 
accommodation, noise reduction, and stability are met by the feedback 
controller. Robustness is a control objective for both controllers. Conversely, SOFFT 
provides a methodology in which the feedforward control law is designed with one cost 
function while the feedback control law is designed with a different cost function. The 
two controllers are combined into the SOFFT control structure (see Figure 
3 above and Figure 4 on page 23) so as to cooperate with each other and 
achieve all the control objectives. 
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H. INTEGRATED FEEDFORWARD AND FEEDBACK (IFF) 

CONTROL DESIGN 


In general, a control system contains both a feedforward and feedback control law. 
We will define the feedforward control law to be the part which depends only on the 
command variables. The feedback control law is the remaining part which depends on the 
plant output with or without random noises and disturbances. For linear syst ems , these 
definitions uniquely specify the feedforward and feedback control laws. For nonlinear 
systems, some ambiguity may remain depending on the particular system. 

With notable exceptions [6], [8], [13], [14], the feedforward control design problem has 
received little attention while the feedback control problem has been studied extensively. 
This asymmetry has resulted in design techniques and methodologies which attempt to 
meet all the control objectives using only feedback control laws. This approach places an 
unnecessary burden on the feedback law which must try to achieve conflicting objectives. 

By an integrated feedforward and feedback (IFF) control, we mean a control law in 
which the feedforward and feedback are designed so as to cooperate with each other. It 
is to be hoped that an IFF control law does not place conflicts between feedforward and 
feedback. 

1. Linearization 

In the SOFFT approach, we often model the plant by a linear system of the form 


Xk+l — &X 4* Fx U xk "t" W x k "t" dx 


( 6 ) 


Vxk = C x Xk + v x k + b x 


( 7 ) 
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where x* is the rc x -component state vector, u x * the n ux -component control vector, y x * the 
rc yx -component measurement or feedback vector; w x k and u x k are Gaussian white noise 
sequences representing system noise and measurement errors, respectively; d x is a constant 
vector which may represent a constant disturbance or the constant term resulting from 
linearizing a nonlinear system about a given operating point. Similarly, b x represents the 
bias error in the measurement coupled with the constant term from the linearization of 
a nonlinear sensor. We may sometimes assume that the actual values of b x and d x are 
selected according to a Normal or Gaussian distribution function. 

In most applications, the plant and sensors are nonlinear systems. The motivation for 
using a linear formulation as shown in (6) and (7) comes from the linearization of nonlinear 
systems. 

The class of linearizable systems plays an important role in atmospheric flight vehicle 
dynamics. For example, the system parameters such as the coefficients of lift and drag are 
defined and determined from given flight condition data; i.e., a linearization about a given 
operating point. The kinematics can be expressed analytically but can also be linearized. 
Thus, the class of linearizable systems is large and contains important application systems. 

Consider the nonlinear system 

dx(t) , 

= i(f) = u x (t), w x {t)) (8) 

Vx{t) = g(x(t), P x (t)) ( 9 ) 

where x(t) is the n x -component state vector, u x (t) is the n uz -component control vector, 
w x (t) is a second-order Gaussian white noise process with zero mean and is assumed to 
be independent of the Gaussian initial condition x(0). The functions / and g are, in 
general, nonlinear functions of their arguments. / describes the plant, actuators, random 
disturbances and sensor dynamics, while g models a nonlinear sensor or feedback vector, 
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y x (t). The measurement noise i/ x (t) is also assumed to be a Gaussian white noise process 
with zero mean, statistically independent of w x (t) and x(0). 

Assume that the system is linearizable or that / and g have Taylor series expansions 

of the form 


d / ^ / 

x(f) = /( 2 0, U X 0,0) + (x O ,U xO ,0) [x(t) - Xo] + (xo,UiO,0) [Ui(t) — U z o] 

+ (x o ,u xO ,0) \w x {t) - 0] + ol(t) (10) 

y x (t) = g{x o,0) + (x o ,0) (x(i) - x 0 ] + ^ (x o ,0) P x (t) + o 2 y {t) (11) 

where xo and u x o are arbitrarily selected state and control vectors respectively, o\[t ) and 
o 2 (t) are the remainder terms of the series and are of 2 nd order in the perturbed variables. 
(xo,ti 2 o) determine a particular operating point or condition. It is important to note that 
/(x o ,«xo,0) need not be zero in this formulation. Thus, both steady-state conditions 
(/(xo,u x o,0) = 0) and unsteady conditions can be accommodated in this formulation. 

In a given problem, some of the components of (xo, «o) may remain constant. Fur- 
thermore, not all the components of (xo,uo) may affect the plant dynamics (/(x, u,u>)) 
in a nonlinear manner. The remaining components of the operating condition (x 0 ,u x0 ) 
which have significant variations for a given problem are included as components in the 
parameter vector, p. Thus, the parameter vector has n p components where 

n p <n x + n ux (12) 

Usually, n p would be much smaller than the right-hand-side (RHS) of (12), as only few 
components of interest are varying for a particular problem. 

The parameter vector p represents the operating point of the system. When designing 
a variable-gain control law for a nonlinear system (see Sections IV and V), the parameter 
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vector, p, plays an important role. Various plant models axe selected to cover a range of 
operating conditions. These plant models correspond to appropriate values of the param- 
eter vector, p. Thus, a clear understanding of the relationship between p and the physical 
variables (xo ? «o) is necessary. 

Now collecting terms and rewriting (10), (11), we get 

x = A x (p) x + B x {p) u x + w x + d x (p) (13) 

Vx - C x (p)x + i/ x +b x (p) (14) 

where the independent variable t has been dropped for notational convenience. The re- 
maining terms can be found by equating terms with (10), (11). 


d f 

A x (p) = T (x o ,Uzo,0) 

(15) 

B x (p) = duT (x o ,u zO ,0) 

(16) 

°x(p) = dxT ( x 0 > 0) 

(17) 

W x [t) = (xo,«xO,0) w x (t) + ol(t) 

(18) 

Vx{t) = di/T (xo,0) u x (t) + oj(0 

(19) 

dx[p) = /( x O,«x0,0) - \A x (p) x 0 + B x {p ) u l0 ] 

(20) 

b x {p) = g{x 0 ,0) - C x (p ) x Q 

(21) 
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Finally, the system in (13), (14) can be described by discrete system [10] when u x (t) 
remains constant from one sample to the next while p remains constant over the sampling 
period, 


Zfc+I = 4>x{p ) Xfc + r x (p) u xk + w xk + d x (p ) (22) 

Vxk = C x (p) x k + u xk + b x (p) (23) 

where <t> x (p) is the state transition matrix and T x (p) the control effectiveness matrix as 
given by [10] . 

If the operating point parameter vector p varies over time, then (22), (23) represent 
a time-varying linear system. When the operating point remains near a given point, the 
system in (22), (23) becomes a time-invariant linear system such as (6), (7). Both cases 
will be considered in the following sections. 

It is important to note that in this formulation, x k , y xk and u xk represent the total 
state, sensor measurements and control commands; they are not the more commonly seen 
perturbed variables. Also note that the plant noise w x contains two effects: random 
disturbances such as wind gusts, electronic or other noises, etc., and second order nonlinear 
effects which come in due to the linearization of the system. 

2. Time-Invariant SOFFT Control 

When the operating point parameter vector p remains near a given operating condi- 
tion, the linearized system given by (22), (23) becomes time-invariant and can be described 
by (6) and (7). In this section, we will summarize the SOFFT approach for linear time- 
invariant (LTI) systems. The details of the development can be found in the Phase I 
report*. 

*Halyo, N., “A Stochastic Optimal Feedforward and Feedback Control Methodology for 
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In the SOFFT approach, we separate the feedforward and feedback control objectives 
into two groups. Then we optimize each group of objectives to obtain each controller. This 
also implies separating the design of the feedforward and feedback control systems as is 
described below. 

The overall control objective is to enable the plant (an aircraft or other 
vehicle or process), to closely track the input commands at every instant in the 
presence of disturbances, despite uncertainties about the plant, within digital 
implementation constraints and other system limitations. The standard regulator 
problem is obtained when the input commands are zero, so that a feedforward controller 
is not needed. The feedback control law stabilizes the plant and maintains the system 
state near zero. The set point regulation problem is obtained when the input commands 
are constant. The most general case where the plant state tracks the input commands 
at every instant is considered in this work. The commands may take the plant from one 
operating point to another, constantly maneuvering in a transient rather than steady state 
condition. 

In the SOFFT approach, the feedforward control has the objective of tracking 
the input commands and producing the desired transient response to com- 
mands (overshoots, damping, etc.) assuming no measurement noises nor plant 
disturbances. On the other hand, the feedback control law has the objectives 
of measurement and plant noise (disturbance) suppression and stability within 
considerations that uncertainties and variations in the plant, actuator and 
sensor subsystems will be present. 

Thus, the feedback gains should be high enough for disturbance suppression consistent 
with low noise feedback and robustness with respect to unmodeled dynamics or parameter 
uncertainty. On the other hand, the feedforward control is not concerned with noise 

Superagility”, ICS FR-689102, Information & Control Systems, Incorporated, 28 Research 
Drive, Hampton, VA, 1989. 
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feedback problems as it is completely implemented on a digital control computer and 
generates a “perfect trajectory” which will be described later. 

Consider the plant given by (6) and (7). Now suppose that we want to achieve 
a trajectory which we denote by {y xk ,k > 0}. Let {u* k ,k > 0} be the control sequence 
which produces the desired trajectory when no random noise (plant or measurement noise) 
is present in the system. Accordingly, we must have 

x k+l = <f> x x* k + T x u* xk + d x (24) 


Vxk — x k 


(25) 


3. Feedback Control Design 

The actual plant output trajectory { y x k,k > 0} will be perturbed by plant distur- 
bances and measurement errors. So that the actual trajectory will differ from the desired 
trajectory in (24), (25). In the SOFFT design methodology, we want to design a feedback 
control law which will keep the actual plant output trajectory y xk close to the desired 
trajectory {y* xk ,k > 0} at all times and despite plant uncertainties. 

Thus, define the tracking error variables 


— %k %k 

(26) 

U xk = U x k — U* xk 

(27) 

Vxk ~ Vxk Vxk 

(28) 


Using (6), (7), (24) and (25), the tracking error variables are seen to satisfy 
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xic+l = 4>x Zk + r* U x k + w xk 


(29) 


Vxk = C x X k + Uxk ( 30 ) 

Now suppose that the desired output y* xk is available at the k th sampling instant t 
Since the actual output y x k is also available from sensor readings, the output tracking error 
y x k is available for use in the feedback control law. 

If u x k can be chosen (designed) appropriately to achieve plant noise or disturbance 
suppression, sensor noise attenuation and stability under plant uncertainties, then the 
tracking error will be small. The input commands, the command model variables or the 
desired *-trajectories are not directly involved in the design of the feedback control law. 
Note that the feedforward control has not entered into the feedback control design. What 
is of interest is to keep the error variables small. This is quite different than the explicit 
model following approach described in the Introduction where the feedback and feedfor- 
ward control laws are designed simultaneously to minimize a single objective function given 
in (4). 

In order to achieve the feedback objectives, the designer may decide to use instanta- 
neous output feedback, control rate commands, reduced-order dynamic compensation, as 
well as integral error feedback. Selecting the feedback control structure which builds the 
control design model is one of the most important parts of the design process. In general, 
the feedback control structure will be of the form 


ujfc-l-l — fife "h At r r Vxk H" W u k 


(31) 


Ck+ 1 = 4>c Ck + 4>cy Vxk + <frcu “fc + At T cx V x k + A tT c V c k + ™ck ( 32 ) 
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Ik+i = Ik + A((H y y,k + Hu it) + w ik 


(33) 


U z k — Cr "t" -^r (34) 

where u* is the n r -component vector of control variables which use a rate command struc- 
ture (i.e., not all controls need have this structure), v x k is the ra uz -component vector 
containing the rate commands for the applicable components while the remaining compo- 
nents are the same as in the original control vector u x jt, c* is the n c -component dynamic 
compensator state vector, v c k is the n vc -component compensator control (design) vector, 
Ik is the n/-component integrator state and At is the sampling period. In (33), H y and 
H u define the variables fed into the integral error feedback subsystem. In general, we will 
use H with a subscript to define tracking error variables and command variables; we will 
use C with a subscript to denote the feedback variables; i.e., the output or feedback vector 
(see (38), (39), (42)). The number of commands is riff. When integral error feedback is 
used, the number of integrators nr equals n# ; otherwise (integral error is not fed back) 
m = 0 and (33) vanishes. The vector sequences {u> u fc}» {tUcfc} and {w/jt} are assumed 
to be Gaussian white noise random processes which can be used as design parameters to 
achieve various objectives. 

When the designer selects to command the rate of change of all the control variables 
in u x k, then 


r r — C r — I , D r — 0 , n r — n uz . (35) 

If the designer does not select a control rate structure, then n r and u* vanish and D r is 
the identity matrix. When the designer selects to command the rate of some control com- 
ponents while commanding the actual position value of the remaining control components 
r r ,C r and D r are correspondingly determined. 
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The dynamic compensator can be forced to be a reduced-order estimator of some 
unmeasured state variables by appropriate selection of the cost matrices Q,R, in the 
objective function. Alternately, the dynamic compensator may be used for other purposes 
such as signal shaping and stability margin. 

Finally, the integrated tracking error can be fed back to achieve a type-1 response for 
the closed-loop system. The tracking error is modeled as H y y x k, and will be discussed 
further later. 

The control feedback structure given in (31) - (34) forms the feedback design model. 
The design model is the augmented system given by 


f%k+l\ / 4x ^xCr 0 0 \ / X/f \ 

«fc+i _ I o I 0 0 | Ujfc | 

Cfc+l I 4cyC x 4cu 4c 0 I Cfc 

\ik+iJ Va tH v C x A tH u 0 iJ \i k J 

/ T x D r 0 \ w x k y 

AiT r 0 | / v x k\ , w u k 

AtTcx AtT c I \v c k ) u>ck 

VO 0 / V t ujk J 

f$xk\ c x 0 0 0 \ /Xfc\ ( Uxk \ 

yuk I _ 0 C u 0 0 Uk U uk I 

Vck I 0 0 C c 0 Ck Vc k I 

\yik J Vo 0 0 cj \i k j \vikJ 


(36) 

(37) 


Xk + 1 — 4 Xk + r vk + w k 


(38) 


Yk = C Xk + Vzk 


(39) 
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The feedback control law is of the form 


v k = -KY k = -KCX k - Ku xk 


(42) 


The feedback gain matrix, K, can be partitioned so as to be compatible with the 
partition of Y k in (40) . 


TV- / Kxx Kx\l Kxc Kxl 

\K CX K cu K cc Ka 

The general formulation described above is highly versatile. You can select practically 
any feedback structure to achieve almost any control objective. Building and optimizing 
the feedback design model is made quite easy by a software product (ACET) developed in 
conjunction with SOFFT methodology. 

The following feedback cost function can be used to optimize the parameters (i.e., the 
gain matrix K ) of the feedback control system. 

J = Jim ■^^E^{xZQX k +2X’[Mv i +vlRq) (44) 

w— oo 2{N + 1) 

Minimizing this cost function requires the use of the Stochastic Output Feedback 
algorithm. 

4. Feedforward Control Design 

In the SOFFT methodology, the objective of the feedforward control law is to track the 
input commands and produce the desired transient response (e.g., overshoots or critical 
damping) when no random disturbances or noises are present. For example, if a fast 
maneuver is commanded, we would like for the control to move as fast as allowed by the 
plant’s physical constraints to achieve the maneuver. 

A high-gain feedback control law has disadvantages such as insufficient sensor noise 
attenuation as well as lack of robustness to plant variations and nonlinearities. Thus, the 
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feedback control law must be selected with care not to result in too high a loop gain. 
On the other hand, the feedforward control law has no such constraint. It is completely 
implemented on a digital computer and except for round-off error, it has no random noise 
or plant nonlinearities to contend with. Furthermore, if a fast and large input command 
is given, it is suppose to produce a fast and large control to produce the response. Thus, 
the feedforward control law is supposed to be a high-gain system, but one which does not 
have the associated disadvantages of high-gain feedback. 

Consider a command model of the form 

Zk+i = <f>zZk + T z u zk (45) 

Vzk = H z z k (46) 

where the command model state z k has ^-components which will usually be different 
(smaller) than the plant model state dimension n x . The command model input, u z , is 
assumed to be a Gaussian white noise sequence with zero mean independent of the initial 
condition z<> Thus, z k and y zk are random processes. 

Two somewhat different interpretations of the command model are possible. In one 
case, the command input u zk is the actual input to the compensated system, such as the 
pilot stick input to the aircraft. In this case, the command model output y zk represents 
the desired output of some plant variables. This case will be used to design a flight control 
system for the F/A-18 aircraft in a subsequent section. 

In the other case, u zk does not represent a physical variable which is measured. The 
command model simply represents a class of responses one of which will be commanded. 
For example, y zk may represent a flight path generated by a flight management computer 
according to certain rules. Accordingly, we can devise a command model {(f> z ,T z ,H z ) 
which is statistically representative of this class of paths. Then u zk can be computed from 
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the given path and the command model selected. In both applications, the design of the 
feedforward control law follows the same procedure. 

We will assume that at time kAt , only the current and past values of the commands 
are known for use in the feedforward controller; i.e., {z*, Zk-i, • • • , zo} and u z k are known 
at the k tfe sampling instant. This constraint is necessary to design a realizable or causal 
control law which can be implemented in real-time. If the commands were known for all 
future values of k, then the feedforward problem would become deterministic. Optimal 
tracking techniques are available for such deterministic cases [15]. However, the result of 
the optimization is an open-loop control sequence rather than a feedforward control law. In 
special cases, it is possible to obtain a control law by making appropriate approximations. 
The CGT [13], [14], [15], belongs to this latter class, as it obtains the feedforward control 
gains by assuming that u z k is a constant vector for all k. However, then it uses this 
feedforward law in cases where u z k varies. This and other approximations in the CGT 
largely account for its unpredictable tracking performance. The SOFFT approach 
avoids making such approximations by casting the problem as a stochastic one 
in which the constraint of causality can be naturally incorporated. 

Recall that we must find the desired control and output trajectories u xk and y xk which 
satisfy (24) and (25) in order to separate the feedback design problem (defined by (26) - 
(30)) from the feedforward control problem. Since the feedforward objective is to track 
the commands, particularly during transient maneuvers, we want to minimize the tracking 
error at every instant. Let 


e k — H y y* x k — B z Zk , k> 0 (47) 

Thus, we would like to maintain e k at very low values for all k > 0. If e k vanishes 
for k > 0, then we have the “perfect tracking” case. Note that this applies only to the 
feedforward since plant and sensor noises will result in less than perfect tracking in the 
actual system. 
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We select a quadratic cost function in e k and include further terms to accommodate 
any objectives other than tracking. 


r = „ lim wwm E £ {ii^tui + *i T q\ *; + <i -;»} 

N^oo 2 (A + 1J " 

Thus, we would like to minimize the cost function J* subject to the constraints 


(48) 


x k + l ~ < t > x x k + rx u xk + 


(49) 


Zk+1 = <t>zZk + T z u z k (50) 

Since u z k is a Gaussian white noise sequence, (48) - (50) pose a discrete LQG problem 
with an unstable and uncontrollable part. The solution to this problem is given in [5], [6]. 
The resulting control law is 


u* xk = -K; x k - K z z k - K u u zk + u* (51) 

where u* is a constant vector depending linearly on d x and b x . As we will use an incre- 
mental implementation for the feedforward controller, the actual values of u x ,d x and b x 
will not be needed. 

Thus, we can determine the desired control trajectory, u* k from (51). Then, using 
(24), (25) we can determine the desired output trajectory, y xk , iteratively. Since we now 
know u* xk and y xk , we can feed these values into the feedback control system to determine 
the feedback control u xk and obtain the total control command to be sent to the plant. 
Thus both feedback and feedforward can be designed separately, each optimizing its own 
set of objectives. 
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5. Algorithm for Feedforward Gains 

To obtain the optimal feedforward gains K *, K z and K u , the necessary conditions for 


optimality may be solved. From [5], these equations can be found to be 

K = [rl p;, r, + p;, (52) 

PL = - r, k lfpu\*, - r , K\ + 1 vj x; k‘ x + h? h, + <?; (53) 

^ = [rf PL r« + -Ri‘] _1 rJ p;, ( 54 ) 

PL = [<L - T, K\ T P^ 4 ,, - Hi H, (55) 

x* = [rj p; x r 1 + «;|-'rj p- x r, (56) 

where 

Hx = Hy C x (57) 


The first two equations (52) and (53) are seen to be the standard full-state linear 
quadratic regulator (LQR) necessary conditions when the state and control penalty ma- 
trices Q and R for the problem are set by 

Q = Hi H x + Q\ , R = Rl (58) 

where Q\ and R\ correspond to (48). A standard LQR algorithm can be used to compute 
K* and P* s . It is important to note that most algorithms will generate a critical point 
rather than the absolute minimum. 
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The necessary conditions for the remaining gains K z and K u are seen to be linear 
equations. P* z is the solution of the discrete Lyapunov equation given in (55). K z and 
K u can be easily computed from (54) and (56), respectively. 


Perfect Tracking 

We define perfect tracking as the case where H y y* k = H z z* for all k and all com- 
mands, u z k . When Q{ and R J vanish, the optimal solution will often produce perfect 
tracking. However, when perfect tracking is not possible, there may continue to exist an 
optimal solution which will approximately track the commands. 

When the matrix [H x r x ] is non-singular, perfect tracking is possible and the feedfor- 
ward gains can be obtained with greater ease. In this case, we have 


K* z = [H x T x ]~ l H x <f> x (59) 

K z = -[H x T x \- l H z <f> z (60) 


K u = -[H x T x ]- l H z T z (61) 

Substitution of these gains into (52) - (56) shows that the necessary conditions for 
optimality are satisfied. While these gains produce the feedforward control necessary for 
perfect tracking, they may produce high gain values and their use should be examined on a 
case by case basis. Thus, when H x T x is nearly singular, the perfect tracking capability of 
the system may be limited due to control authority limitations. In such cases, a different 
tracking variable may be selected. 
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6. The SOFFT Control Structure 

The control command sent to the plant is the sum of two components, one generated 
by the feedforward controller, u* xk , the other by the feedback controller, u x *.. The feedfor- 
ward control u xk is determined by (24), (25), (50) and (51). The feedback control u x * is 
determined by (31) - (34) and (42). An incremental implementation is recommended for 
the actual digital implementation of the SOFFT control law. 

Combining both feedforward and feedback control laws results in a system structure 
depicted by the block diagram in Figure 4. This diagram shows the feedforward control 
law in greater detail. Considering only the structure, of the SOFFT control law, rather 
than the methodology, we note that the SOFFT structure has 

1) feedforward dynamic compensation beyond the command model, 

2) two feedforward links to the feedback loop: u xk and y xk 



NOISE 


Figure 3. SOFFT integrated feedforward and feedback (IFF) control structure 

It is also important to note that these two feedforward links (u^,!/**) are highly 
correlated. As can be seen from Figures 3 and 4, the correlation is introduced by the plant 
and sensor models used in the feedforward control system. When these models perfectly 
match the actual plant and sensors, then the actual and desired outputs y x k and y xk match 
at every instant and no feedback correction u x k is needed. When the feedforward models 
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Figure 4. SOFFT Control Law Structure 










do not match the actual plant and sensors, whether due to random noises or modeling 
errors, the outputs y x k and y* zk do not match and the feedback control system produces a 
corrective action, u x k- Thus, the feedforward and feedback systems are integrated so that 
they cooperate in trying to achieve all the control objectives. 

It is important to note that neither the error feedback structure shown in Figure 1 
nor the explicit model follower structure shown in Figure 2 match the SOFFT structure. 
In comparison to the Command Generator Tracker (CGT), note that SOFFT uses a dy- 
namic compensator while the corresponding CGT feedforward control law consists of a 
gain matrix. Also note that the CGT and the model follower structures do not feed 
directly. 

The SOFFT control law does not require full-state feedback. It accommodates tiny 
sensor configuration by using stochastic output feedback. Further note that the number 
of controls and the number of outputs or feedback variables do not have to match, but are 
arbitrary. 

Finally, the SOFFT structure appears to be a combination of the model follower and 
error feedback structures in that it feeds forward both y* k and u zk . 
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III. COMMAND TRACKING PERFORMANCE 


Most current measures of robustness treat the sensitivity of the feedback control law. 
For example, the classical phase and gain margins provide an excellent measure of stability 
robustness for single-input single-output (SISO) systems. The loop gain and phase may 
vary by the respective margins before the closed-loop system becomes unstable. For multi- 
input multi-output (MIMO) systems, similar stability margin criteria have been developed 
using <7-plots or the singular values of the return difference matrix, cr-plots also provide a 
measure of the error in a MIMO error feedback control system. 

While the analysis of the feedback loops provides invaluable information, it analyzes 
only one part of the total system. When the system has a feedforward control law reaching 
the feedback loops at more than one point, it is desirable to have a measure of the error 
in tracking which can handle the total system. 

In this section, we will develop a measure of command tracking performance (CTP) 
or the ability of the compensated system to track input commands. We call these analysis 
tools ^-plots. 

First we develop a mathematical model of the complete SOFFT control law including 
the feedforward and feedback control laws in their general form. Here we will allow for 
the plant and sensor models used in designing the feedforward control law to be different 
than the ones used in designing the feedback control law, both of which may be different 
than the actual plant and sensor model. We will denote the feedforward models with the 
superscript “ * ” while the actual models will not have this superscript. 

There are a large number of state vector selections which can be reasonably selected 
to describe the system. All of these are valid. Here, we have selected the following form 
of the feedforward and feedback systems compensating the plant. 
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1 *k+i = 4> z zk + r* u zk 

x* k+l = [r s -T: k:\xi -T* s K z z k 

(62) 

FEEDFORWARD 

- r x K u u xk + {dl + T* x u* x ) 

(63) 

CONTROLLER 

u\ k = ~K x\ -K z z k - K u u zk + u* t 

(64) 


\ y* xk = C:x* k + D* s ul k + C* z z k + C:u zk + bl 

(65) 

PLANT | 

j %k+ i = fix x k + V x u xk V x u xk ■+• w xk + d x 

(66) 

(PADS) j 

\ Vxk = c x x k + D x u xk + u xk + b x 

(67) 

( 

t Ufc +1 = U k + At r r v xk + w uk 

Ck+l — 4>cC k + 4>CU U k + <t> C y y xk + At T c V ck 

(68) 

FEEDBACK 

1 + AtT cx v xk + w ck 

(69) 

CONTROLLER / 

j ^ ^ 

h+ i =h + A t[Hy y xk + H u u k ] + w Ik 

(70) 


Uxk = C r u k + D r v xk 

(71) 

CONTROL | 

i 

1 v k = ( v , xk ) = —K x [y xk - y**] - K u u k 
! \ v ckJ 


1 

— K c c k — Ki f k 

(72) 

SYSTEM j 

1 H y y xk = H y C x x k + Hy u xk + H y b x 

(73) 

OUTPUT | 



TRACKING j 
ERROR / 

Hy Vxk = Hy [ yxk 2/xfc] 

(74) 


where the gain matrices K X ,K U ,K C and Kj are partitions of the feedback gain matrix K 
in (42) and (43) . 

v — i v v is v \ — ( H xx K xu K xc K x j \ 

K-(K X K u K, K,)-y K " ^ K " K 'J 

Note that although D x and D* have been neglected (assumed null) in the development 
of the SOFFT feedforward control law, we are including them in this formulation to obtain 
a more general analysis tool. Also note that some notational changes have been made in 
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comparison to our Phase I final reports, such as v uk and T cu have been renamed v xk and 
r cx , respectively. 


After some manipulation the system equations (62) - (74) can be put in state space 
form as follows. 


z k + 1 > 


x k + 1 
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C* = C* 2 +D x K z 
C* = C* + D x K* 
C* = C* U + D x K u 


( z k\ 

x k 

*k 
Uk 

C It 

\h) 


D x — D x D x 
b x = b x - bi 


(76) 


(77) 


where the constant forcing terms related to d x ,b x ,d x , etc. and random noises have been 
neglected as they will not impact the desired transfer functions. 

Using the state-space model above, frequency response from the input command u z k to 
the outputs H y y xk and H x y* xk can be computed with any selection of the parameters, with 
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the feedback loops open or closed. In these computations, it is assumed that D x D r = 0. 
The outputs are obtained using 


Zk \ 
x k 

*k ( -D X K U 

s* V ft 

Cjfe 

Ik) 


c; = c;-d;k z 
c; = c:~ D* K* 
c: = c:-d;k u 


1 . Z-Plots 

Let T(u>) be the transfer function matrix from the pilot command input u z (u > ) to the 
plant output or feedback vector, y z (w); i.e., 

y*(w) = T(u)u z {(jj) . ( 80 ) 

Let the desired or commanded plant outputs, y*(w), be given by 

y*(w) = T*(w) u z {u) , (81) 

When we have more sensors (plant outputs) than commands, the corresponding trans- 
fer functions from u z (u>) to the commanded output are H y T(u ) and H y T*(u). We as- 
sume that there are as many pilot inputs (n u J as there are commanded plant outputs; 
i.e., H y T{u) is n Ui x rt Ux . We will further assume that H y T*(u > ) is invertible, a.e. 

We define the Z-matrices as 


= ( 82 ) 




( V*k\ _ ( -D.K, -D X K‘ C x D X C, 0 0 ) 
W*J V < 5 ; 6 ; o o o o) 
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z 0 (w) =[ff v r(w)| [s,r («)]-* 


(83) 


To understand the meaning of the Z-matrices, first note that 

Zj(u>) = Zo{u) = I , when H y T(u) = H y T*(u) . (84) 

Thus, when the actual response equals the desired response, the Z-matrices axe unity. 
Now define the u,(w),tt*(a/),y(u;),y*(w) as follows, 

(85) 

( 86 ) 


u*(u>) = Z/(w)u z (w) 


y{u) = ZoM y*(u;) 


These definitions lead to the following interpretations when T(u>) ^ T* (cj); i.e., when 
off-nominal conditions are being analyzed. 


u 


* 

Z 



»*(“) 


y‘H 


Figure 5. Zi-Plot interpretation 

The system response to a pilot input varies form the ideal response in off-nominal 
conditions to achieve the desired response y*(u>), the pilot constantly varies his input 
commands u z (u>). If the pilot has to vary his inputs drastically, the system is not robust. 
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Conversely, if only small adjustments in the pilot input are needed, the system’s tracking 
performance is robust. Note how this is related to pilot workload. 

Figure 5 shows the effect of Z/(cj) in terms of block diagrams. It is seen that if the 
pilot compensates for the off-nominal condition by including Zf 1 ^) as a transfer function, 
he will produce perfect tracking. Thus, both the magnitude and the frequency content of 
ZfHu) is important as a measure of the amount of work the pilot must do to compensate 
for off-nominal conditions. 







Figure 6. Zq- Plot interpretation 


While Z / (u/) is a measure of the input sensitivity, Zo(u>) is a measure of the output 
sensitivity. If a pilot input u*(w) produces the output [H y y*(w)] = y*(w) with the nominal 
aircraft, then the same pilot input will result in the response [H y y x (w) ] = y(u>) when 
the aircraft is in off-nominal conditions. Thus, Zo[<J) measures the sensitivity of the 
commanded plant outputs to off-nominal conditions when the pilot makes no effort to 
adapt to these conditions. Figure 6 shows the effect of Zo(&) in terms of block diagrams. 
When the control system (feedforward and feedback) does not vary, Zo (cj) depends only 
on the feedback loop parameters (plant and feedback controller). On the other hand, Z/(u>) 
depends both on the feedforward and feedback loop parameters. We will refer to Z/(u>) 
and Zo{<jj) as the ^-matrices and refer to the analysis plots obtained from them as Z-plots, 
for convenience. We think that the ^-matrices and the transfer function T-matrices, T(w), 
provide highly useful information for the analysis and design of MIMO control systems. 
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Now define the following measures 


Z ma x{u) = || Z((j) || 2 

(87) 

Z min {u) = \\Z~\u)\\i' 

(88) 

= 1 + || Z(u) - diag {£(<*/) }|| 2 

(88) 


Recall that ideally the Z-matrices (and their inverses) should equal the identity matrix. 
This can be seen from Figures 5 and 6 and from (84). Thus, for a robust MIMO system, 
the diagonal elements of the Z-matrix should be close to 1, while the off-diagonal elements 
should be small. As defined above, Z max (u>) and Z m i n (u) bound the ^-matrix above and 
below, whereas Z c (u> ) provides a measure of the magnitude of the off-diagonal elements. 
Thus, Z c (u>) isolates the magnitude of cross-talk; e.g., the impact of command 1 on output 
2. Note that for a SISO system, Z c (u>) vanishes as there are no off-diagonal elements. 

These measures are conservative. They display the worst case scenario for all the 
loops. It is often important to look at the individual elements of the ^-matrices and 
the T-matrices. This is particularly necessary when the system loop characteristics are 
different; e.g., when the desired bandwidth for two loops is significantly different. Thus 
(u/) and Ta (w) are also important. 

2. F/A-18 Actuator Approximation Design 

To illustrate the techniques developed a longitudinal flight control system design for 
the F/A-18 was obtained at a single flight condition. The intent was to test the methodolo- 
gies developed for approximation of actuators (or more generally the plant model) during 
the feedforward controller design, and also to test the usefulness of the sensitivity measures 
developed. 
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These considerations led us to a MIMO SOFFT design. The flight condition used 
was level flight at 20° angle-of-attack. The controls were the throttle and stabilator. The 
sensors used were a pitch rate gyro, an angle-of-attack vane and the thrust level. Although 
three sensors are fed back, the pilot can command only two variables independently since 
there are only two independent controls. The commanded plant outputs were selected 
to be: 1) q - H |a (i.e., a pitching command), 2) ST (i.e., a thrust level command). The 
command model for the pitching command was chosen to be the short-period mode ap- 
proximation. For the second command, a \ at order model with 2 sec. time constant was 
used. 


I 




a) Pitching Response 


h) Thrust Response 


Figure 7. Desired response (command model) 


Using the design approach developed, a digital feedback control law with a PIF* 
structure was designed at a 25 Hz sampling rate. The feedback design included actuator 
♦The PIF structure consists of a proportional feedback portion, an integrated error feed- 
back portion and a filter portion which weighs the control rate. Thus, the PIF structure 
can filter out high frequencies, produce type 1 behavior and adjust time constants. It has 
been successfully used in various control laws including [14], [17], [25]. 
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models for both the thrust and stabilator. On the other hand, in the feedforward design, 
the stabilator actuator was neglected from the plant model. The throttle-to-thrust model 
was kept since this has a long 2 sec. time constant. 

The desired response to pitching and thrust commands is shown in Figure 7. All plots 
are shown in db for amplitudes and degrees for phase. Two feedback designs were made. 
Feedback (FB) Design 1 was the first design made. Using a “perfect tracking” feedforward 
design (without stabilator actuator), a SOFFT controller was obtained. Figure 8 shows 
the Z-plots for this design. From Figure 8a, we see that the pitching and thrust response 
will be within 1 or 2 db of the ideal or desired response for any pilot command until 40 
rad/sec. However, some cross-talk is apparent from Figure 8b. 

Figures 8c and 8d show a higher level sensitivity. If the pilot were to try to compensate 
his input commands so as to achieve the ideal response, significant adjustment would be 
needed in the higher frequency range above 1 rad/sec. 

After some trial-and-error, we obtained Feedback (FB) Design 2. The Z-pIots for this 
design are shown in Figure 9. While the ||i?o ||2 has been reduced modestly, the “loop 
cross-talk” has been reduced impressively for both Zj and Zq. The cross-talk from the 
pitching command u Zl to the thrust response is responsible for the difference. 

From Figure 10, it is seen that the maximum level has been reduced by 40 db, from 
-20 db to -60 db near 2 rad/sec. Actually, depending on the design goals, both designs 
display characteristics. However, the Z-plots seem to show specifically what’s good and 
what’s not as good. We axe working on some small extensions of the concepts. 

The individual elements of the ^-matrices are not shown here. Our experience so far 
with Z-plots and T-plots have been rewarding. These measures seem to provide important 
information about the tracking performance of the system. The sensitivity of the command 
tracking performance of a SOFFT design will be analyzed using Z-plots in section VI. 
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a) FB Design 1 JJ FB Design 2 

Figure 10. Frequency response from pitching command (u*,) to thrust response 
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IV. VARIABLE-GAIN SOFFT DESIGN METHODOLOGY 

In the preceding sections, we have dealt with linear time-invariant systems. Of course, 
the time-invariant or single-model SOFFT control law tends to be very robust. So that it 
can control a nonlinear plant over a wide range of operating points even though the design 
model corresponds to a single operating point. 

As the actual operating point moves away from the design condition, the assumptions 
used in the design methodology become less valid. And although the performance of the 
single-model SOFFT control law may still be acceptable, it is clear that a design based 
on the current operating point can have higher performance. Furthermore, by adapting to 
the current operating point, we can extend the operating regime of the plant indefinitely 
within the physical constraints of the plant. Thus, extending the operating range 
and having greater performance throughout the operating range are the major 
motivating factors for the Variable-Gain SOFFT design methodology. 

For aircraft, the operating point is the flight condition of interest. To produce high 
levels of agility or maneuverability, it is necessary to move from one flight condition to 
another with great ease and high flight path accuracy. When the flight conditions are 
relatively close, a robust control system can provide satisfactory performance within a 
small flight regime. However, when the flight regime extends beyond certain limits, a 
constant-gain control system no matter how robust cannot maintain the necessary levels 
of agility. Or conversely, a control system with a variable-gain structure can provide greater 
agility when appropriately designed. 

In this section, we will develop the design methodology for Variable-Gain SOFFT 
control laws. Given the successful applications [23], [24], [27] of the Variable-Gain Output 
Feedback design methodology [17], [18], we expect that Variable-Gain SOFFT will also 
provide a successful method for extending the flight regime with improved flying qualities. 
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1. Problem Formulation 


In section II. 1, the nonlinear plant and sensor models in (8) and (9) were linearized 
and discretized resulting in the discrete system given by (22) and (23). 

When the operating point varies, it is necessary to linearize the plant about a new 
operating point (x 0 , u 0 , 0) or p at every iteration. As long as the sampling interval, say 
At, is selected sufficiently small, the operating point parameter vector p remains near its 
initial value over the interval. Thus, the linearization about the new plant, say Pk, is valid, 
and the original nonlinear system can now be described by the plant and measurement 
model 


Zfc+1 = <t>i(pk) *k + r z (pjfe) u x fc + w x k + d x (pk) 


(90) 


VxI k = c x (p k ) Xk + u xk + b x (p k ) 


(91) 


To design a control law for this system, consider the following class of time-invariant 
systems. 


x(p,k + 1) = <j> x {p) x(p, k) + r z (p) u x (p, k) + w x (p, k) + d x (p ) (92) 

Vx{p, *) = C x (p) x(p, k) + u s (p, k) + b x (p) (93) 

where p is a parameter vector of n p components which represents the particular flight 
condition or operating point. For a fixed p, (92) and (93) represent a linear time- invariant 
system. The parameter vector p can be chosen by the designer in the way that best suits 
the particular problem under investigation. In particular, p may be a nonlinear function of 
the state vector or the measurement vector; or it may represent an unmeasured parameter 
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of the system. In all cases, the vector p must be either measured or estimated using a filter 
or other parameter estimation technique. 

The system matrices r(p), etc. can be arbitrary functions of p; however, they 

must be known functions of p. For example, we can define these matrices by specifying 
them at several critical flight conditions, and then defining them by linear interpolation in 
between these flight conditions. Of course, the matrices can also be defined analytically 
when such expressions are available. Finally, note that the variables x(p,jfc), y x {p,k), 
represent total quantities rather than perturbations of the physical variables. This 
is achieved by appropriate interpretation of the terms d x (p ) and b x (p). 

» In the Variable-Gain Control formulation, we shall allow the command model to vary 

with the operating point parameter vector p. 

i . i 

z{p, k + 1) = 4> z (p) z(p, k ) + r 2 (p) u z (p , k) (94) 

1 

e(p, k) = Hy(p) y x (p, k) - H z (p ) z(p, k) (95) 

The ability to vary the command model with the parameter vector provides an im- 
portant flexibility in the control system. First, it allows the designer to command different 

| 

variables (or different linear combinations of the variables) as the operating point or flight 
condition enters significantly different modes of operation; e.g., you can command different 

2 

variables during normal flight than you would during stall or post-stall operation. 

The second flexibility is that the commanded handling qualities can change with the 

i 

. . s 

operating point parameter vector p. Thus, the desired response to a given command can 

be different according to the value of the parameter vector. Also recall that the parameter 

vector can be selected arbitrarily. Thus one component of p may be chosen as a mode 

1 selector to represent fast or sluggish handling. Then according to the value set by the 

pilot, the aircraft would produce a fast or slower response. There are clearly a wide 
variety of ingenious ways of defining the parameter vector to produce desirable effects. 
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The objective is for H y (p,k) y x (p, k) to track the commands H z (p,k) z(p,k) at every 
sampling instant. Alternately, the objective is to minimize the tracking error e{p,k) at 
every sampling instant. 

Following the SOFFT approach for single model time-invariant systems, suppose that 
the control trajectory u* x (p,k ) produces the desired feedforward state trajectory x*[p,k). 
Then 

z‘(p,k + 1 ) = Mp) **(»>.*) + r.(j + <*«&>) < 96 ) 

Now consider the “local” cost function J*(p): 


J*(p) = - lim — Y' Ei\\e*{p,k)\\l + x*{p,k) T Q\{p) x*{p,k) 

2 iV—oo N + 1 “ t 
k = 0 

+ u* x {p,k) T R{(p)u* x {p,k ) | ( 97 ) 

where e*(p, k) is defined in accordance with (47). 

When the feedforward objective includes only command tracking, the weighting ma- 
trices Q\ and R\ can be set to zero, so that only the tracking error is weighted'. In other 
cases, the feedforward objectives may include additional goals such as keeping the control 
co mman ds u* x at relatively low values. These can be included using the additional terms Q* 
and R \ However, the added objectives will be obtained at the expense of some additional 

tracking error. 

The cost function J*(p ) describes the local objective at the operating point or flight 
condition corresponding to the vector p. Since the real objective is to design a control 
system for the complete operating range or the flight regime, a global cost function can be 
defined by a linear combination of the local costs over the region of interest. 


flmod 


r=Y, A-Hp') 


1=1 


(98) 
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where p‘ denotes a particular value of the parameter vector p and determines a particular 
operating point in the operating range, /» is a weighting coefficient which can be selected 
by the designer to specify the relative importance of the particular local operating point 
p * and n mo d is the number of operating points included in the cost function. Note that 
nmod should be chosen so that the complete operating range is represented satisfactorily. 

Following the variable-gain output feedback methodology, let us constrain the form of 
the gain matrices as shown below. 


Kif) = K> + E pU) Kj (99) 

3 = 1 

Tip 

K z (p) = K z0 + Y^p(j)K zj (100) 

3 = 1 

Tip 

Ku (p) = K u0 + ^2 p{j) K uj (101) 

3 = 1 

where p(j) is the ] th component of the parameter vector p, K* } , K ZJ and K UJ are constant 
gain matrices. The control law is of the form 

u *x{Pi k ) = -K*(p) X*(p,k) - K z (p) z(p,k) - K u (p)u z (p,k) + u*(p) (102) 

Thus, the problem is to minimize the global cost J* in (98) subject to the constraints 
of (96), (94) and (102). An algorithm to obtain the gain matrices in (99) - (101) will be 
described in section IV.3. 

Perfect Tracking 

To achieve perfect tracking, the cost function is selected so as to weight only the 
tracking error e*(p, k). Thus, the matrices and in the local cost J*(p) are set to 
zero. In this case, the local cost can be minimized subject to the local plant model at 
each operating point corresponding to p. Thus, it is not necessary to limit consideration 
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to a finite number of operating points, p*, as in the cost (98). In this case, the complete 
operating range can be included. The solution is the same as the time-invariant case except 
that now the plant parameters are functions of p. 

K*(p) = [H x {p)T x (p)\~ l H x {p) <j> x {p) (103) 

K z (p) = —\H x {p) rx(p )]- 1 H z (p) Mp) (104) 

K u (p) = -\H x (p ) r x (p )] _1 H z (p) T z {p) (105) 

It is important to note that the optimal gain matrices for the perfect tracking case 
are not of the form shown in (99) - (101). Of course this is not significant as long as the 
gain matrices can be computed with ease in real time. 

Since the plant model matrices are known a priori, the only computational load in 
obtaining the gain matrices is in inverting the matrices [H x (p) r x (p)]. The dimension of 
this matrix is n yx by n UI . Thus, it is necessary to invert a matrix of dimension as large 
as the number of controls (n uz ) or the number of commands (n yi ). When the number 
of controls is small such as 1, 2 or 3, it is possible to invert the matrices analytically. 
For higher order matrices, inversion routines are available. However, as the dimension 
increases, the computational load increases as the cube of the dimension and may reach 
undesirable levels. Also, for high order matrices, the accuracy of the inverse may be called 
into question. However, for up to third order, the inversion can be accomplished with ease. 

2. Feedforward Integrators 

For single-model designs, it is often possible to achieve perfect tracking in the feedfor- 
ward control. As before, by perfect tracking we mean that the desired output trajectory 
matches the command model output at every sampling instant; i.e., 
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e* (p, k) = H y (p) y* (p, fc) - JJ, (p) z(p, A) = 0 , k > 0 (106) 

In other cases, perfect tracking may not be achieved. For example, in a single-model 
design in which > 0 (iEJ is positive definite), perfect tracking will not be obtained even 
though the tracking error may be small. In the Variable-Gain feedforward formulation 
given in the previous section, perfect tracking will not be achieved in the general case. For 
non-minimum phase systems, it may be desirable not to use the perfect tracking solution. 
Finally, in other cases, the plant’s physical constraints such as position and rate limits may 
make it more desirable not to seek perfect tracking in the feedforward system as well. 

Thus, in general H y yl(p,k ) will not track the commanded H z z(p,k ) perfectly at 
every k and p. As long as the tracking error is small, this may not be significant. However, 
in some cases, it is desirable to have no error in steady-state for constant commands. For 
this purpose, it is possible to introduce integrators into the feedforward control law. 

Thus, let 


I* (p> A: + 1) = /* (p, A:) + At e* (p, k) 


(107) 


<(P> k ) = -KiP)^(P^)-^z{p)z{p,k)-K u (p)u z (p,k)-K* I {p)r(p,k) + u* x (p) (108) 

Now, we introduce an integral penalty term into the objective or cost function. Now 
consider the local cost function 

J ' = \ FTT ^ + Q\{p) x’(p,k) 

k~0 ' 

+ u* x (p,k) T r; (p) u* (p, k) + r (p, k) T Qj I* (p, k) | (109) 

while the global cost remains as given by (98). 
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The result of minimizing the new cost function will be a type-1 behavior in the feedfor- 
ward system. In the transient, some tracking error will be present, but the error will vanish 
in steady-state. The algorithm for gain computation is given in the following section. 


3. Algorithm for Gain Computation 

In developing an algorithm for computing the feedforward variable-gain matrices, we 
will include the case for feedforward integrators. To exclude the integrators, simply neglect 
the corresponding terms 

First consider the augmented co mm and model 




( 110 ) 


( 112 ) 


Z(p,k + 1) = I j>z(p) Z(p,k) + r z (p) w z (p,k) (111) 

where w z (p,k) is the white noise sequence whose covariance equals that of u.(p, fc). 

Mp) = (*i p] T i p) ) , r z (p) = (°) 

Now, consider the augmented state vector X*(p,k). 

( . 

X‘(p,k) = \r(p,k)\ ( 113 ) 

\z(p,k) 

After some manipulations, it follows that 


X*(p,k+ 1) = <j>{p) X*{p,k) + T(p) u*(p, k) + Mp) w z (p,k) 


( 114 ) 


Y'(p,k) = X*(p,k) 

where the augmented matrices can be built as given below 


(115) 
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4>{p) = 


(116) 


( <Mp) 0 0 \ fr m {p)\ 

A<tf x (p) I —AtHz(p) 1 , r(p) = 0 

V 0 0 <f> z (p) J \ 0 J 


C = I , D = 0 , T w (p) = ( 0 ) , H z (p) = {H,(p) 0) (117) 

\Tz(p)J 

Now we express the local cost function in terms of the augmented state. 

(H?(p) H x (p) + Q\(p) 0 -ffj(p) Hz(p)\ 

Q‘ (p) = o Qi(p) 0 (118) 

V -ffj(p) H x (p) 0 Hl(p)Hz(p) j 


J 'tP) = \ j™ jA~L ^ E{x‘(P'k) r Q‘(p)X‘(p,k)+ul(p,k) T R; u;(p,*)| (119) 

fc =0 ^ ' 

Umod 

J‘ = £ ( 120 ) 

t=l 

Observing (114) - (120), note that we have embedded the stochastic optimal feedfor- 
ward control design problem into the Variable-Gain Output Feedback problem treated in 
[17] and [18]. 

Thus, the algorithm for the current problem is obtained by building the augmented 
matrices in (116) - (118) for each of the models considered; i.e., for each value p‘, 
i = 1, 2, ..., n mo d. Then by using the Variable-Gain Output Feedback algorithm in [18] 
pp. 28-29, we obtain the optimal feedforward gain. 


u x(P» fc ) = ~K*{p)X*[p,k) (121) 

K-(p) = (K- X (p) KUp ) K,(p) K4p)) = K- 0 + £ pU) K- (122) 

;'=i 
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Our experience with this algorithm to date indicates a fast rate of convergence, in 
fact, much faster than usual variable-gain output feedback problems. We think this is due 
to the fact that C = /; i.e. , the feedforward uses a full-state feedback structure. 

4. Digital Implementation 

The implementation of the SOFFT feedforward control law is intrinsically different 
than that of feedback control law. In the feedback control law, only those equations 
involved in computing the control commands, u xk , are implemented. In particular, the 
plant state and sensor output vectors are not computed; the actual plant and sensors 
produce the outputs which are input to the feedback control system. 

In the SOFFT feedforward control law, both the control commands, u* fe , and the 
desired trajectory, y* xk , are computed in the implementation. Usually, this involves an 
approximate plant and sensor model. Of course, the command model must also be imple- 
mented. 

We have developed an incremental implementation obtained by differencing and some 
approximations. The usual advantages of incremental implementations are that trim values 
are not needed, that integrator wind-up problems do not occur and limiting is easily 
handled. 

Feedforward Control Implementation 

Let pk be a parameter vector value at the sampling instant t k . 


&Ugk — U ? k - Uzk-l 

(123) 

Azfc +1 = 4>z{pk) + r 2 (pA;) A U z k 

(124) 

zjfc+1 - Zk + Azfc+1 

(125) 
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a <* = -Kin) a *; - k;w MKM vU-x - KM **-,] 

Kz{pk) AZfc -^u(pfc) AUzk 


(126) 


Axfc + 1 = <t>* x {pk) Axfe + r;(p fc ) A u* xk (127) 

Ay* jk+i = C*(p k ) Ax£ +1 (128) 

y*fe+i = vlk + &ylk+i (129) 

4 = (130) 


Thus, the feedforward control law produces the vectors u* xk and y* k +i which are now 
input to the feedback control law to obtain the complete or combined control law. The 
gains K*(pk),K}(pk),K z (pk ) and K u (pk) are computed on line from (99) - (101) or from 
(103) - (105) for perfect tracking. Note that K}(pk) = 0 for perfect tracking; otherwise, 
it has the form shown in (122) or (99) - (101). The plant and command model matrices 
are also computed on-line using the interpolation algorithm described below. Of course, 
other interpolations may also be used. 

It can be shown that the perfect tracking control law coupled with the above incre- 
mental implementation will produce the desired y* k trajectory such that 

Hy{Pk) Ay lJt+ i = H z (pk) Az^i , k > 0 (131) 

When H y and H z cure independent of p, this results in the perfect tracking trajectory 


Hy y x k+i — H z Zk+ 1 > k > 0 


(132) 
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provided that the initial conditions match; i.e., 

H y y * 0 = H z z 0 (133) 

This produces a type-1 property for steady-state tracking as k — * oo. On the other 
hand, the same type-1 property does not hold for the optimal tracking option. A constant 
command may not produce the equality shown in (132) in steady-state (or as k — ► oo). In 
many cases, this type-1 property may not be necessary as long as the transient behavior 
is satisfactory. To ensure steady-state tracking with zero error, it is necessary to include 
integrators in the feedforward control law. 

Interpolation of Feedforward Models 

As mentioned earlier, the feedforward control law generates the *-trajectories using 
a model of the plant and sensors. In the variable-gain design case, these models are 
specified at a finite number of operating points sufficient to represent the operating range. 
In actual operation, the operating point parameter pk will naturally have values different 
than the specified models. To accommodate the complete operating range, it is necessary to 
interpolate in between these models. The following interpolation algorithm was developed 
for this purpose. 

Let p 3 , for j = 1, 2, • • • , M, represent the } tfl model parameter vector. Let p be the 
current value of the parameter vector. Define the metric pj{p) as follows: 

Pj{p) = WP-VW > j = l,2,..M (134) 

where || • || is the 2-norm for vectors. Of course, any norm could be used for this purpose. 
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Algorithm 


1. Find the n smallest values of Pj{p),j = 1,2, If Pj{p) = 0 for some j, then 

<£*(p) = stop. It may be desirable to also check that the selected n models are 
not too close to each other; e.g., ||p* — p*\\ > e. 

2. Reorder the models so that the selected models are the first n models, 

3. Compute 

1 




Note: 


v m _ , 

Mp) 


(136) 


4. Compute the interpolated matrices <£*(p),r*,C*(p) using the formula 


«"-£ S'- 


XX > 


K< = Kip') 


(137) 


Note: When n = 2,n p = 1, this algorithm produces the usual linear interpolation. 
Feedback Control Implementation 

The feedback control uses the sensor measurements t/ x * and the feedforward variables 

Vlk’<k, etc - 


At Jxk — Vxk Vxk—1 


(138) 


V x k — [J — K xu ^pif) C U T r ]v x k —1 Kxx{Pk) [Ay z fc Ay x jg] ^xc{Pk] C c (pfc) Acfc 

- A tK x i{pk) Ci{p k )[Hy{pk ) Vxk-1 ~ Hy{p k ) yjjfc -1 + H u {pk) Uk-1 ] (139) 
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Ac * +1 — ( 0 C (pfc) - At T c {pk) K cc (pk) C c {pk)) Ac k 




I 

i 


! 


+ (0cy(pfc) - At T c (pk) K cx (pk)) [Ay x k - Ay* fc ] 

- {At) 2 r c (p fc ) K cI {p k ) Ci(p k ) \Hy[pk) Vxk-i - H y [pk) y* xk -i + H u [p k ) ufc_i] 
+ (0c«(pfc) — At r c (p fc ) -K" c U Cp*) Cu{pkj) At r r v xk — i 


+ Afr cz (p fc )(u lfc -V xk -i) (140) 

Cfe + l = Ck + Ack+l (141) 

Ufc = u k _i + At T r Vxk-i (142) 

«I k = C r UJ k + Dr Vxk + U*jfc (143) 


The control commands u xk are now sent to the plant. When p* is constant, the 
implementation reduces to the single-model case. 


| 


1 


| 
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V. F/A-18 SIMULATION 


A general dynamic flight simulation module was developed to test the new control de- 
sign methodology proposed in this work. This simulation module incorporates the specific 
aerodynamic properties of the F / A-18 aircraft, although any other aircraft can also be sim- 
ulated by inserting the particular aerodynamic and propulsive modules for that aircraft. 
The non-linear F/A-18 simulation was also integrated into ACET for various platforms. 
This combination of ACET and non-linear aircraft simulation produces a very powerful 
tool to design flight control systems. 

1. Overview of the Basic Components 

Figure 11 identifies the basic components of the F/A-18 simulation and the interac- 
tions among them. For each integration period, At, the aircraft simulation receives the 
command inputs from the flight control system (FCS), computes the total forces and mo- 
ments on the aircraft by using engine, actuator and the aerodynamic models of the aircraft, 
obtains the position and speed components using the aircraft kinematics and produces sen- 
sor outputs to be used by FCS. Each dynamical model in the simulation can be improved 
or simply changed by varying the model dynamics while keeping the structure of the simu- 
lation intact. The actual algorithm of the simulation also involves various data exchanges 
among the basic components. The simulation is intended for arbitrary flight maneuvers in 
relatively short periods of time, so the body axes components axe selected to be computed 
and the Earth-fixed axes are an inertial frame of reference. The simulation also contains 
a vectored thrust capability. The control actuator models have linear l at -order dynamics 
with non-linearities such as rate and position limiters. The kinematics equations for a rigid 
airplane are solved using a non-linear Runge-Kutta-Fehlberg method for the single step 
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Figure 11. Flight simulation diagram 
























integration. Since F/A-18 can be simulated in different mission conditions, arbitrary cen- 
ter of gravity locations in the aircraft can be selected. Steady components of wind velocity 
vector are also included. An ideal sensor model is employed for sensor parameters. The 
F/A-18 simulation was integrated into ACET for VAX computers on the VMS platform 
and for IBM-compatible computers with DOS. Combining the model building, control de- 
sign, analysis, linear simulation, matrix operations, and plotting capabilities of ACET and 
the non-linear F/A-18 simulation, the integrated product enables the researcher to increase 
his efficiency in designing flight control systems greatly. Detailed important features of the 
F/A-18 simulation components are: 

Flight Control System (FCS) Inputs 

The FCS module provides the simulation section the necessary pilot input variables 
such as the control surface actuator commands, engine throttle and vectored thrust com- 
mands. The FCS inputs to the simulation as elements of the command input vector, U 
and their references are: 

17(1) = Stabilator (positive trailing edge down (+TED), deg) 

U{ 2) = Aileron (+TED, deg) 

17(3) = Rudder (+TE left, deg) 

U[4) = Throttle angle (deg) 

17(5) = Pitch vane command for vectored thrust (VT) 

U (6) = Yaw vane command for VT 

Any combination of the command input variables can be selected with the simulation 
setting the other input parameters to constant values. The command input variables are 
computed by FCS for each integration interval, A t. 

Actuator Dynamics 

Actuator dynamics model includes a linear l st -order dynamics for the actuator sur- 
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faces and non-linearities consisting of position and rate limiters. The actuator commands 
from FCS are used in a configuration module to obtain the positions of the actuator sur- 
faces. A list of the control actuator surfaces and their references are: 

DHTR = Right stabilator (positive trailing edge down, (+TED)) 

DHTL = Left stabilator (+TED) 

DAR = Right aileron (+TED) 

DAL = Left aileron (+TED) 

DRR = Right rudder (+TE left) 

DRL = Left rudder (+TE left) 

DTFR = Right trailing edge flap (-(-TED) 

DTFL = Left trailing edge flap (+TED) 

DLFR = Right leading edge flap (+TED) 

DLFL = Left leading edge flap (+TED) 

Speed brakes (0 - 60 deg) (set to 0) 

Landing gear (0-1 down) (set to 0) 

Due to asymmetric position limiters, the stabilators, the ailerons and the all flap surfaces 
may have asymmetric responses. 

Engine Dynamics 

The non-linear engine model for the F/A-18 aircraft was incorporated in the simula- 
tion. The engine model is based on data received from DFRF and its implementation in 
AGCB at LaRC. The engine commands from the FCS include the throttle command which 
determines the various conditions of the engine such as the flight idle, military power, af- 
terburner on, etc. A vectored thrust capability using pitch and yaw thrust commands also 
exists. 

The engine dynamics are modeled in an engine module which requires angle-of-attack, 
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mach, altitude and attitude rates and produces total engine forces and moments. The 
iteration time for the engine is 0.032 sec and it is synchronized with the F/A-18 simulation 
by using additional iterations if necessary. The engine module accounts for various losses 
such as ram, inlet/aft body, windmill drags and thrust vectoring losses to compute the net 
thrust of the engine. The gyroscopic effects are also included in computations. 

Static and dynamic tests have been made with the engine module at various throttle 
level, mach number and other conditions. The rate limit, gyroscopic and other non-linear 
effects as well as phugoid and short period modes due to attitude changes have been 
observed in the engine output and are found reasonable. 

Aerodynamic Model 

The aerodynamic model of the F / A-18 aircraft in the simulation computes the aerody- 
namic forces and moments using modules supplied by Langley researchers. The model uti- 
lizes the control surface actuator positions from the actuator models and the current body 
axes components from the kinematics modules. The model includes a quasi-static-elastic 
approximation where the flexible modes are actively modeled to compute the aerodynamic 
coefficients. The force and moment coefficients and the lift and pitch derivatives are first 
computed about the aerodynamics reference center (a wind tunnel reference), and then 
translated into the body axes with respect to center of gravity (CG) location. The model 
also includes a rolling moment increment for a > 40° as well as other effects for high a. 
Various mission conditions can be employed changing CG location and mass distribution. 
The simulation employed a default configuration representing an early estimate of F/A-18 
HARV. The forces and moments obtained from the aerodynamic and the engine models 
produce the total forces and moments to be used in the kinematics model. 


55 


Kinematics 

The kinematics formulation for the simulation was developed using a six-degree-of- 
freedom formulation. The equations of motion were derived for a rigid airplane and a 
quasi-static-elastic approximation was used to compute aerodynamic forces and moments. 
The motion of the aircraft was computed in the body axes, rather than the stability axes 
since the simulation was developed for various flight maneuvers rather than a particular 
flight condition. For the moment equations, symmetry about the x-z plane was assumed. 
The simulation was primarily intended for maneuvers occurring in relatively short periods 
of time, rather than long periods, such as cruise and navigation. Consequently, a flat, 
non-rotating Earth model was assumed, so that the Earth-fixed axes are an inertial frame 
of reference and the mass and inertia distribution was assumed to be constant during the 
simulation not dependent on factors, such as fuel load, weapon adjustments, etc. For 
each simulation, various mission conditions can be selected with different mass and inertia 
distributions and center of gravity locations. Gravity is assumed to be constant with the 
default value selected at sea level. It is also constant over the airplane volume acting as a 
force at the center of gravity with no moments produced. 

The kinematics equations used in the simulation are given below: 


u = —q xv + r v — g sin 0 - I F x 

771 


(144) 


v = pxv — r u + g cos & sin <f> + — F v 


(145) 


w — q u ~ pv + g cos 0 cos H F z 

m 


(146) 


P = K , 


( hh-Jl !k ) tr + K iP , +{j-)L+(-^)N 


(147) 


4 = ( J -i^F) r p + (ifi) (r 2 - p 2 ) 4- (f ) M 

ly 1 y 


(148) 


56 



r = Kx 


(It ~ Izly ~ Igx\ v- , r - 1 A r 

(— Yt ^-)pq-K 2 qr + ——L+- r N 

lxlz •**■** -«z 


4> = p + sin <p tan 9 q + cos (j> tan 9 r 


9 = cos <p q — sin <f> r 


ij) = sin 4> sec 9 q + cos (f> sec 9 r 


(149) 

(150) 

(151) 

(152) 


x = (cos 9 cos ip) u + (sin <p sin 9 cos tp — cos <p sin <p ) v 
+ (cos <p sin 9 cos ip + sin <p sin ip) w 


(153) 


y = (cos 9 sin ip)u + (sin <p sin 9 sin ip + cos <p cos ip) v 
+ (cos <p sin 9 sin ip — sin <p cos ip) w 


(154) 


z — — sin 9u + sin <p cos 9 v + cos <p cos 9 w 


(155) 


where 


K i = 


T 2 

-(h) 


Ki = lh (I* - r v+ «] 


In these equations, the body axes origin is at the center of gravity and the axes are 
oriented such that +x is out the nose, +y is out the right wing and +z is out the belly. 
The airplane is ass um ed to be symmetric about the x — z plane, therefore, I xy and 'y* 
products of inertia are zero and are not included in the kinematics equations. 
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The kinematic equations are solved using a non-linear single-step integration algorithm 
employing a seventh-eight (7,8) order Runge-Kutta-Fehlberg method and a solution for 
all kinematic parameters are obtained at each simulation step size, At. The integration 
module employs variable internal step sizes depending on the variability of the solutions 
and adapts to existing singularities and invalid inputs. The simulation options include 
selecting relative and absolute accuracies. The algorithm was found to be performing well 
with good numerical stability. 


Sensor Outputs 

The F/ A-18 simulation was developed to give the option of selecting a number of sensor 
outputs necessary for the flight control system. An ideal sensor model is employed and any 
combination of sensor outputs can be fed back to the control system. The available sensor 
parameters include all the computed kinematic parameters, body angles, body mounted 
accelerometers, etc. A complete list of available sensor outputs are given below: 


y(i),...,y(i2) 
Y{ 13) 
Y{ M) 

y( 15 ) 

y(i 6 ) 

y(i 7 ) 

y(is) 

y(i9) 

y(2o) 

y(2i) 

y(22) 


= X(l) X(12) From kinematic equations 

= Vt, True air speed (ft/sec) 

= a, Angle-of-attack (deg) 

= 0, Sideslip angle (deg) 

= Fbx, Body mounted accelerometer in the x-axis (ft/sec 2 ) 
= FbYz Body mounted accelerometer in the y-axis (ft/sec 2 ) 
= Fbn, Normal body mounted accelerometer (ft/sec 2 ) 

= h, Sink rate (ft/sec) 

= 7, Flight path angle (deg) 

= p, Barometric pressure (lbs/ft 2 ) 

= p, Air density (slugs/ft 3 ) 
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2. Simulation Examples 

F/A-18 aircraft simulation had been tested in various static and dynamic conditions. 
Various modes including phugoid, short period and lateral modes had been analyzed; static 
tests were compared with the trim conditions given by Langley researchers at various mach 
numbers and attitudes. Dynamical tests had also been made by using non-trim initial 
conditions for the aircraft dynamics. 

Figure 12 shows the longitudinal/vertical variables where longitudinal variables have 
perturbed initial conditions. The control inputs are set for a trim condition at a = 18.8° 
and mach number = 0.24. The plots indicate a short period mode with a period of 8 or 9 
sec., and a phugoid mode with a period of about 40 sec. Both of these natural frequencies 
are in general agreement with the theoretical values determined from the linearized models. 
Figure 13 shows the variables in a simulation of the same flight condition where lateral 
variables have non-trim initial conditions. The natural frequencies of these modes are also 
in general agreement with the theoretical values. 


ACET™ ACET™ 




(a) Airspeed (b) z&-axis speed 

Figure 12. Longitudinal/ vertical variables with longitudinal variables perturbed 
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ACET™ 



ACET™ 



(c) Pitch rate 


ACET™ 


(e) Negative altitude 


(d) Pitch angle 



Figure 12. Longitudinal/vertical variables with longitudinal variables perturbed 
(continued) 
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SIMULATION PLOT 


SIMULATION PLOT 


(c) Lateral speed 


(d) Yaw rate 


Figure 13. Simulation variables with lateral variables perturbed 
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ACET™ 



(e) Yaw angle 


ACET™ 



SIMULATION PLOT 
(f) Pitch rate 


ACET™ 



(g) Pitch angle 


Figure 13. Simulation variables with lateral variables perturbed (continued) 


B 
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Figure 14 shows the effects of the nonlinear engine dynamics in the simulation. The 
aircraft is essentially in trim with level wings as initial condition where the control inputs 
are set at a = 18.8° and mach number = 0.24. A large throttle pulse to maximum throttle 
and back to its trim position is applied as shown in Figure 14(a). The engine gross thrust, 
the x-axis thrust and pitching moment are shown in Figures 14(b) - 14(e). Note that the 
gross thrust oscillates with the phugoid mode, but the x-axis thrust seems more stable. 
Also note that the gross thrust is shown for one engine, while the net x-axis thrust includes 
the effect of both engines. Rate limit and other nonlinearities are seen in Figure 14(c) which 
zooms on the leading edge of the pulse. 
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(a) Throttle position input 



ENGINE SIMULATION 
(b) Gross thrust 


Figure 14. Engine variables with a maximum throttle pulse input 
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(c) Gross thrust (zoom to see detail) 


(d) Net x-axis thrust 
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Figure 14. Engine variables with a maximum throttle pulse input (continued) 
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u (rr/s) 


In Figure 15 the simulation plots show that the throttle command excites all the 
modes of the aircraft dynamics. The phugoid mode is clear in Figure 15, the short period 
mode rides over the phugoid in Figures 15(d) and 15(e). The gyroscopic effects seem to 
produce a significant response in the lateral variables as seen in Figures 15(f) - 15(i). The 
combined response of engine and airframe seem reasonable. 



SIMULATION PLOT SIMULATION PLOT 

(a) Airspeed (b) Pitch angle 


Figure 15. Simulation variables for the same maximum engine throttle pulse input 
as in Figure 14 
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(g) Yaw rate 


(h) Lateral speed 


ACET™ 



(i) Roll angle 


Figure 15. Simulation variables for the same maximum engine throttle pulse input 
as in Figure 14 (continued) 
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VI. A VARIABLE-GAIN SOFFT FLIGHT CONTROL SYSTEM 

FOR THE F/A-18 


The SOFFT control design methodology developed in the preceding sections can be 
applied to any nonlinear system which can be linearized and discretized as described in 
section II. This class of linearizable systems is quite large and includes very diverse applica- 
tions from flight controls to robotics, from crystal growth to power transmission networks. 

Here we will illustrate the use of the SOFFT design methodology by applying it to 
a flight control system (FCS) design for the F/A-18 aircraft. The objective is to enable 
the aircraft to provide high performance with high flying/handling qualities throughout 
its flight envelope. Alternately, we want to extend the flight regime by improving the 
aircraft’s flying qualities until the physical limitations of the aircraft are reached, such as 
control authority, rate limits, etc. 

1. Longitudinal FCS Design 

To demonstrate the methodology, a variable-gain SOFFT control law was developed 
to design an angle-of-attack (a) command system for the F/A-18 aircraft. A single-model 
a-command system was also designed in the earlier stages of methodology development. 
This will be used to compare the effect of using variable gain. 

The aircraft longitudinal/ vertical dynamics was modeled by a 4 tfc -order linearized 
state-space model of the form 

x = A(p) x + B(p) 6 a + w(p) + d{p) (156) 

x T — {u' a q 0) (157) 
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where u' is the forward speed along the z body-axis in units of 10 ft/sec, a is the angle- 
of-attack in degrees, q is the pitch rate in deg/sec and 0 is the pitch angle in degrees. The 
control is the stabilator (<5 a ) or horizontal tail position in degrees. 

Feedback 

Both the feedback and feedforward control laws used variable-gain. For this purpose, 
3 linearized models corresponding to 3 flight conditions in level flight, with level wings at 
three different angles-of-attack and airspeed were used. The operating point parameter 
vector p was chosen to be the angle-of-attack (a). The 3 flight conditions were selected at 
5°, 20° and 40°. Thus, 

p = a , p 1 = 5 , p 2 = 20 , p 3 = 40 . (158) 

For the feedback design model, the plant model was augmented by an actuator system 
for the stabilator. The stabilator actuator dynamics were assumed to be the same for the 
various flight conditions. 

5 a = —30 6 a + 30 u x (159) 

Using a sampling period (At) of .04 sec (i.e., a sampling frequency of 25 Hz), the 
augmented system was discretized. 

The output or feedback vector y x was selected as 

»* = (a) + "* +4 * (160) 

Thus, a pitch rate gyro and an a-vane were used as sensors for feedback. The main reason 
for using these sensors is that they are body-mounted, basic measurements. Whereas a 
pitch gyro may not be as reliable in many maneuvers. Of course, since p is needed for 
the adaptive variable gains, a measurement of a is desirable. The only drawback of this 
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Figure 16. Frequency-domain loop analysis of a-command variable-gain feedback 
control law: model 1; perspective: input. 
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Figure 18. Frequency-domain loop analysis of <*-command variable-gain feedback 
control law; model 3; perspective: input. 
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sensor configuration is that the phugoid mode cannot be adequately damped. However, 
since pilots generally produce phugoid damping often without conscious effort, this is not 
a significant setback. 

A Proportional-Integral-Filter (PIF) structure was used for the feedback control law. 
This structure produces a type-1 closed-loop control system as well as allowing control rate 
weighting to achieve the necessary feedback objective of noise suppression. Of course, other 
feedback structures can also be used to achieve the desired feedback objectives. For the 
current problem, the PIF structure adds 2 more states (the control u s and the c*-integrator 
I) to the longitudinal dynamics and actuator. Thus, the feedback design model has a total 
of 7 states. A digital Variable-Gain Output Feedback [17], [18] control law was designed 
using the ACET CACSD tool. 

Feedforward 

In the SOFFT approach, the flying qualities are largely achieved by the feedforward 
control law. This is due to the fact that the SOFFT feedforward law produces the initial 
system response which is highly weighted in pilot ratings. On the other hand, in the error 
feedback structure, the feedback produces both the initial and steady-state responses to a 
given command. 

Whereas the stabilator actuator was used in the feedback design model, it was not 
included in the feedforward design model to reduce the order of the feedforward control 
law. Thus, the 4 tfc -order discretized longitudinal dynamics model was used. It is important 
to note that a short period approximation of 2 nd -order can also be used if further order 
reduction is desired. However, the designer must use iT-pIots to evaluate the consequences 
of such approximations in terms of command tracking performance. 

The command model was selected to produce desirable short period flying qualities. 
Short period handling qualities are rated “satisfactory” in a small region of the complex 
plane in the vicinity of a natural frequency of w n = 3 rad/sec and a damping ratio of = 0.7 


73 



ACET™ FREQUENCY ANALYSIS 



FREQUENCY (RAD/SEC) 

Figure 19a. Effect of pitch rate gyro errors on angle-of-attack (a) for Model 1. 



o 

UJ 

& 


LLl 

tfi 

< 


FREQUENCY (RAD/SEC) 

Figure 19b. Effect of ct- vane errors on angle-of-attack (a) for Model 1. 
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Figure 20b. Effect of a-vane errors on angle-of-attack (a) for Model 2 
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(e.g., see [21] pp. 511-518). Thus, the command model was obtained by discretizing the 
following continuous model at a sampling period of 40 msec. 

<**(*) _ yz{f) _ K w _ 3 7 

u z {s) u*(s) s 2 + w n s + w 2 

where u z is the pilot input command and y z is the desired response of the angle-of-attack. 

Thus, a 2 nd -order co mm and model was used. While this methodology allows the 

command model to vary with the flight condition parameter, p, we kept it constant. 

Using the methods described in section IV, we designed a perfect tracking SOFFT 

feedforward control law. Using the incremental implementation described, we simulated 

the variable-gain SOFFT control law. 

2. Lateral FCS Design 

A single-model SOFFT control law was designed to compensate the lateral dynamics. 
With the variable-gain design methodology having been demonstrated by the longitudi- 
nal/vertical a-command system, the main motivation was to close the loop around the 
lateral dynamics to perform a complete simulation. A roll rate and sideslip command 
system was designed for this lateral control law. 

A 4 tfc -order lateral model was used to design both the feedforward and feedback control 
laws. 



x = A x + B u + w + d (162) 

xT = (/? p r <}>) , u T = (5A SR) (163) 

where ft is the sideslip angle in degrees, p is the roll rate in deg/sec, r is the yaw rate 
in deg/sec, 4> is the roll angle in degrees, 6 A is the aileron position in degrees and SR 
is the rudder position in degrees. Note that the aileron variable 8 A corresponds to the 
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anti-symmetric operation of the right and left ailerons surfaces SAr and SAl, respectively, 
whereas SR corresponds to the symmetric operation of the right and left rudder surfaces, 
SRr and SRl, respectively. Thus, 

SA l = 6 A , 6A r = -SA (164) 


SRi, — SR , SRr = SR 


(165) 


The leading edge and trailing edge flap positions were set by an automatic program 
as a function of Mach number and angle-of-attack. 


Feedback 

While the simulation incorporates actuator dynamics models for the aileron and rud- 
der surfaces, these actuators are very fast. For the design of the lateral control law, the 
actuator dynamics were neglected. For the feedback design model, the flight condition 
corresponding to 20° angle-of-attack was used. 

Using a 0.04 sec sampling period, the plant was discretized with the standard sampled- 
data discretization method. 

The lateral output or feedback vector y x was 


y* = 



+ u x + b x 


(166) 


w 

Thus, a sideslip or /3-vane and lateral rate gyros were used as the sensor configuration. 

A Proportional-Integral-Filter (PIF) structure was used for the feedback control law. 
For the current problem, the PIF structure adds 4 more states to the lateral dynamics 
model. The additional 4 states correspond to the aileron and rudder surface position 
commands, one roll rate (p) integrator and one sideslip (/3) integrator. Thus, the lateral 
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feedback design model has a total of 8 states. The digital Stochastic Output Feedback [l], 
[2] technique was used to design the lateral feedback law. 


Feedforward 

The 4 tfc -order dynamics model used in the feedback design model was also used to 
design the feedforward control law. Thus, the feedforward design also neglected inclusion 
of actuator dynamics for the simplicity of the control law . 

The command model for the lateral feedforward used a 2 nd -order system for the roll 
rate command and a l a *-order model for the sideslip command. For the roll rate, the 
command model is the discretized version of 

p * if! = Mifl = ^ f = .7 , w n = 3. (167) 

u*i(s) u*i(s) s 2 +2$w n s + tv 2 

For a unit step input, this model reaches the e -1 level in less than 0.9 sec and reaches 
the co mm anded level of unity in about 1.5 sec, as it produces a small overshoot. Thus, it 
may be compared to a time constant between .5 sec and .9 sec. Pilot ratings of the open- 
loop roll mode are found satisfactory for time constants T R under 2 sec with corresponding 
values of initial roll acceleration (e.g., see [21] pp. 520-523). Thus, the command model 
seems appropriate to achieve the desired response. For the sideslip command model, a 1 
order model with a time constant of 4 sec was used. Thus, the complete lateral command 
model was of 3 r<i -order, and was obtained by augmenting the roll rate and sideslip models 
described above. 

Using the method described in section IV, a perfect tracking feedforward control law 
was obtained. The incremental digital implementation of the lateral control law was used 
in the simulation runs shown later. 
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TABLE 1 


EIGENVALUE ANALYSIS: MODEL 1 
OPEN-LOOP EQUIVALENT s-PLANE EIGENVALUES 


REAL 

IMAGINARY 

DAMPING 

RATIO 

NATURAL 

FREQUENCY 


0.000000 

0.000000 

0.000000 

0.000000 


CRC 

0.000000 

0.000000 

0.000000 

0.000000 


INTEGRATOR 

- 0.002215 

-0.068608 

0.032264 

0.068644 ] 



- 0.002215 

0.068608 

0.032264 

0.068644 j 

rnU vjrL/lJJ 

- 0.547966 

-1.344528 

0.377412 

1.451903 } 


SHORT 

- 0.547966 

1.344528 

0.377412 

1.451903 J 


PERIOD 

-30.000000 

0.000000 

1.000000 

30.000000 | 

STAB. ACTUATOR 


TABLE 2 

EIGENVALUE ANALYSIS: MODEL 1 
CLOSED-LOOP EQUIVALENT s-PLANE EIGENVALUES 


REAL 

IMAGINARY 

DAMPING 

RATIO 

NATURAL 

FREQUENCY 

- 0.004555 

-0.096194 

0.047297 

\ PHUGOID 

- 0.004555 

0.096194 

0.047297 

0.096302 J 

- 1.547369 

-0.449355 

0.960326 

1.611295 | SHORT 

- 1.547369 

0.449355 

0.960326 

1.611295 f PERIOD 

- 5.583268 

-4.599254 

0.771844 

7.233672 

- 5.583268 

4.599254 

0.771844 

7.233672 

-61.775017 

0.000000 

1.000000 

61.775017 
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TABLE 3 


EIGENVALUE ANALYSIS: MODEL 3 
OPEN-LOOP EQUIVALENT s-PLANE EIGENVALUES 


REAL 

IMAGINARY 

DAMPING 

RATIO 

NATURAL 

FREQUENCY 


0.000000 

0.000000 

0.000000 

0.000000 j 

CRC 

0.000000 

0.000000 

0.000000 

0.000000 j 

INTEGRATOR 

- 0.009877 

-0.148653 

0.066297 

0.148981 \ 

[> PHUGOID 

- 0.009877 

0.148653 

0.066297 

0.148981 J 

- 0.251082 

-0.637923 

0.366245 

0.685557 1 

' SHORT 

- 0.251082 

0.637923 

0.366245 

0.685557 J 

PERIOD 

-30.000000 

0.000000 

1.000000 

30.000000 \ 

STAB. ACTUATOR 


TABLE 4 

EIGENVALUE ANALYSIS: MODEL 2 


CLOSED-LOOP EQUIVALENT s-PLANE EIGENVALUES 



! ■ 


REAL 

IMAGINARY 

DAMPING 

RATIO 

NATURAL 

FREQUENCY 

■ ■ - 


- 0.025907 

-0.160355 

0.159491 

0.162434 

J - PHUGOID 

- 0.025907 

0.160355 

0.159491 

0.162434 

- 0.820359 

0.000000 

1.000000 

0.820359 


- 1.342804 

-1.750519 

0.608642 

2.206288 

) SHORT 

- 1.342804 

1.750519 

0.608642 

2.206288 

f PERIOD 

-10.839646 

0.000000 

1.000000 

10.839646 


-41.800201 

0.000000 

1.000000 

41.800201 
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TABLE 5 


EIGENVALUE ANALYSIS: MODEL 3 
OPEN-LOOP EQUIVALENT s-PLANE EIGENVALUES 


REAL 

IMAGINARY 

DAMPING 

RATIO 

NATURAL 

FREQUENCY 

0.000000 

0.000000 

0.000000 

0.000000 

CRC 

0.000000 

0.000000 

0.000000 

0.000000 

INTEGRATOR 

- 0.028034 

- 0.028034 

-0.183087 

0.183087 

0.151357 

0.151357 

0.185221 ) 
0.185221 J 

► PHUGOID 

- 0.279930 

-0.665594 

0.387681 

0.722064 | 

SHORT 

- 0.279930 

0.665594 

0.387681 

0.722064 j 

PERIOD 

-30.000000 

0.000000 

1.000000 

30.000000 j 

STAB. ACTUATOR 


TABLE 0 

EIGENVALUE ANALYSIS: MODEL 3 
CLOSED-LOOP EQUIVALENT s-PLANE EIGENVALUES 


REAL 

IMAGINARY 

DAMPING 

RATIO 

NATURAL 

FREQUENCY 

- 0.067430 

-0.181984 

0.347443 

0.194074 ) 

- 0.067430 

0.181984 

0.347443 

0.194074 } PHUgoid 

- 0.530053 

0.000000 

1.000000 

0.530053 

- 2.460401 

-2.044951 

0.769048 

3.199281 ) SHORT 

- 2.460401 

0.044951 

0.769048 

3.199281 / PERIOD 

- 4.119399 

0.000000 

1.000000 

4.119399 

-36.188185 

0.000000 

1.000000 

36.188185 
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3. Analysis and Simulation 

The Variable-Gain SOFFT control law designed for the longitudinal/vertical dynam- 
ics of the F /A-18 was evaluated in detail. The evaluation of this a-command FCS included 
both analysis and simulation. Note that both the feedforward and feedback control laws 
are direct digital designs whereas the nonlinear aircraft dynamics and actuation systems 
are continuous. So that the closed-loop or, more precisely, the compensated system is 
a sampled-data system. The frequency response analyses are computed directly in the 
z-plane and shown as a function of the frequency u, where z = e _wAt . Thus, no approxi- 
mation back into the w-plane is made in order to analyze the compensated system. Thus, 
the exact sampled-data system is analyzed and then displayed in the format of an analog 
system which is familiar to most designers. 

Eigenvalues 

The eigenvalues of the open-loop plant and the plant with the loop closed by the 
feedback control law are shown in Tables 1, 2 and 3 for each of the models used. These 
eigenvalues determine the stability of the system at the given models. They also show the 

damping ratio of each mode. 

The main objective here is to provide adequate damping for the short period mode 
and the other modes in the feedback compensator. As discussed earlier, while some im- 
provement in the phugoid damping is achieved in each model, sufficient damping is not 
expected from this feedback law as it does not feed back the phugoid variables: pitch and 

airspeed. 

The short period mode has been brought to comfortable levels of damping for a = 
5°, 20° and 40°. 

Feedback Loop Analysis 

The loop gains at stabilator input are shown in Figures 16 - 18 for each model. The 
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Bode plots show that the closed-loop system is in fact a type-1 system and will produce 
zero error in steady-state for constant commands. For model 1 at a = 5°, the Bode plots 
also show a significant drop in the loop gain near .1 rad/sec. This is clearly due to the 
phugoid mode. 

The stability margins can be seen from both the Nyquist and Bode plots. The phase 
margins are greater than 44° for each model. The gain margins are below -10 db for all 
models. 

The c-plots or the return difference are positive except over a narrow bandwidth with 
the lowest value being above -4 db. 


Sensor Noise/Error Suppression 

4 

One of the objectives allocated to the feedback law by the SOFFT approach is to sup- 
press measurement errors. Or at least not to amplify these noises except when necessitated 
by a more important objective. 

To analyze the impact of sensor errors on the aircraft response, we show the frequency 

\ 

response of the angle-of-attack to measurement noises introduced by the pitch rate gyro 
and the a- vane. These plots axe shown in Figures 19 - 21 for a = 5°, 20°, 40°. 

The pitch rate gyro errors are significantly attenuated throughout the spectrum for 

all three models. It is important to note that any bias errors present in the rate gyro will 

i 

have no impact on the angle-of-attack. 

i 

i 

Errors and noises in the a-vane are also suppressed. However, the attenuation comes 
largely in the higher frequencies. Within the command bandwidth which we consider to 
be the low-pass band below 2 rad/sec, a-vane errors are essentially passed through. Thus, 
a-vane noise is suppressed through a low-pass filter. This, of course, is necessary if we axe 
to pass through the largest part of the sensor signal or the true value of the angle-of-attack. 
Thus, a bias error in the a-vane will result in the angle-of-attack being off by the same 
amount. 
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Plant Disturbance Accommodation 


Another important objective allocated to the feedback law by the SOFFT approach 
is to accommodate or suppress plant disturbances. For example from the linearization 
equations in section II (e.g., (10), (11)), we see that when the aircraft moves from one 
flight condition to another, the second order terms may start from zero, but they get 
larger and settle at a non-zero constant value which acts as a plant disturbance. For 
example, the 2 nd order term o\{t) in (10) is zero when the initial condition is (xo,u IO ) 
because o^(t) is a function of the perturbation in x and u. As the perturbations (x — xo) 
and (u z — u z o) get larger, o\(t ) will act as a disturbance. 

Here we have analyzed a plant disturbance acting at the stabilator input such as a 
nonlinearity in the actuator system or the aerodynamic load on the horizontal tail. Figures 
22 - 24 show the effect of such a disturbance on the pitch rate and angle-of-attack responses 
at a = 5°, 20°, 40°. 

Note that significant suppression of the disturbance is being achieved at all models. 
Thus, the impact of most nonlinear effects or other disturbances such as winds is likely 
to be small according to the plots. Note that at a = 5° high frequency artifacts mainly 
outside the command bandwidth will be slightly amplified. This is not present at the other 
a values. A small modification of the feedback law should cure this effect. In general, the 
system appears to accommodate most disturbances. 

Command Tracking Performance 

The preceding considerations such as stability, damping, sensor noise suppression 
and plant disturbance accommodation are the objectives allocated to the feedback control 
law. However, they do not guarantee satisfactory flying qualities without an appropriate 
feedforward design. 

Here we analyze the impact of the feedforward control law on the closed-loop system. 
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Figure 31c. F/A-18 SOFFT control law simulation: a-command = 30°. 
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Thus, we shall look at the response of the angle-of-attack to a command input by the pilot. 

Figure 25a shows the desired frequency response, T*(uj), of the angle-of-attack to a 
pilot input. Thus, this is the frequency response of the command model which has been 
selected to produce flying qualities deemed satisfactory by the pilot. Observing 
we note that it is essentially flat until 1 rad/sec at which point it starts to move down 
at 40 db/decade. We may consider the command bandwidth to extend from 0 rad/sec 
till 2 rad/sec. Thus, we would like to pass through commands within this bandwidth and 
attenuate signals outside. This is why the feedback law can attenuate high-frequency noise 
signals. 

The actual frequency response of the combined feedforward and feedback control laws 
is shown in Figures 25b, c and d for models 1, 2 and 3, respectively. Note the high accuracy 
of these response. In large part, this is due to the fact that we are using a variable-gain 
SOFFT control law rather than a single-model SOFFT control law. Clearly, for small 
pilot inputs, the response of the angle-of-attack will be nearly perfect producing high pilot 
ratings. 

We can assess the command tracking performance mere clearly by looking at the 
system’s Z-pIots, defined in section III. Figures 26 - 28 show the Z-plots for the 
SOFFT control law at each model. Since there is one control, Zj — Zq. Furthermore, 
^max(w) = Z min {uj) - ||Z(w)|| 2 = Z-norm and Z c (u) = 0. Figures 26b, 27b and 28b show 
the weighted ^-norm which weights the Z - norm by the square magnitude of the desired 
response. The intent here is to weigh the command bandwidth more than points outside 
the bandwidth. Note that the Z-nowns for each model are within 2 db inside the command 
bandwidth. The weighted norms show good performance throughout. 

Tracking Performance Sensitivity 

The adaptive nature of the variable-gain SOFFT controller makes it highly robust, 
simply because the control gains are constantly adapting to the current conditions. To 
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assess the robustness of the command tracking performance, we must introduce uncertain- 
ties. 

Now suppose that the operating point parameter vector p used by the control laws 
has a value of 20° while the true value is 30°. Clearly such large discrepancies are not 
likely to occur, even though small errors in the sensor are generally present. However ,it 
illustrates the robustness of the variable-gain SOFFT control law. 

Thus, Figure 29 shows the frequency response, the Z - norm and the weighted Z- 
norm for this situation. Note that the Z - norm is within 3.5 db throughout the command 
bandwidth, but shows higher departures outside the bandwidth. 

Figure 30 shows the scenario when the control law is configured for p = 40° while the 
true value is 30°. Within the command bandwidth, the Z-norm is within 1 db of perfect 
tracking. 


Simulations 

The variable-gain SOFFT a-command control law and the lateral p//?-command sys- 
tem were used to compensate a highly accurate nonlinear simulation of the F/A-18 aircraft. 
The simulation is described in section V. In all the runs shown, the aircraft is initially at 
trim in level flight, with wings level and at 18.8° of angle-of-attack. The throttle (PLAC) 
is set at 60.8° and remains at this level throughout the simulation. 

Figure 31 shows the response to a pilot step input commanding 30° angle-of-attack 
at time 0. The command model produces the desired response for z\. The SOFFT feed- 
forward control generates the desired or *-trajectory y*; y*j is the pitch rate, y* 2 the 
angle-of-attack. Now y* is fed into the feedback loop to form the error signal, y x . Si- 
multaneously, the feedforward generates the command u* shown in Figure 31c which goes 
directly to the stabilator actuator. The feedback control law uses the output error y x to 
generate the feedback portion of the control u x which added to the feedforward command 
u* to produce the total command to the stabilator actuator system. 
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Comparison of the a response (y l2 ) and the command y* 2 shows that they are practi- 
cally the same! The pitch rate response is also very accurate. From the stabilator (DHTL) 
response, note that the actuator has reached its rate limit and the aircraft is pitching up 
at the maximum rate possible for this condition. 

Figure 32 shows the simulated responses to a pilot step input of 40° of angle-of-attack. 
The a response is just as fast; however it has a sharp break point at 39° at 2 sec and levels 
off, reaching 40° at 5 sec. The underdamped phugoid mode is noticeable. 

Figure 33 shows the simulated responses to a pilot step into of 5° at time zero. At 8 
sec, the pilot commands 55° of angle-of-attack. The a response slightly undershoots with 
a breakpoint at 7.5° in under 1 sec. From 5° the angle-of-attack reaches 55° in less than 3 
sec, displaying an excellent overall response. The pitch rate reaches a peak of 45 deg/sec 
in less than 1.5 sec. 

It should be noted that the response is obtained without the benefit of any throttle 
action helping in the pitching moment. Also note that the aircraft physical limit for 
sustained angle-of-attack is about 57° ; so that after 45° the aircraft may be considered in 
the stall region. 

Figure 34 shows the simulated responses for a pilot command of 18.8° of a. However, 
the lateral pilot input commands a 10 deg/sec roll rate at 2 sec for 2 sec, then a -10 deg/sec 
roll rate at 10 sec for 2 sec, while a sideslip of 0.0 deg is commanded throughout. The roll 
rate follows the command model output with accuracy. The angle-of-attack is maintained 
very near its commanded value of 18.8° despite the roll maneuver. Also note the yaw rate 
command and response to produce zero sideslip, and the corresponding drop in the pitch 
angle while a remains constant. 

Figures 35 and 36 show the simulated responses for a constant-gain SOFFT control 
law designed for oc = 20° . While the robust behavior is impressive, the variable-gain 
control performs better. 
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CONCLUSIONS 


The Stochastic Optimal Feedforward and Feedback Technique (SOFFT) is a new 
control system design methodology in which the feedforward and feedback control systems 
are designed separately and then combined so as to cooperate with each other. Traditional 
design techniques optimize a single cost function (which expresses the design objectives) 
to obtain both the feedforward and feedback control laws. In the SOFFT approach, two 
cost functions are defined. The feedforward control law is designed to optimize one cost 
function, the feedback optimizes the other. 

The main feedforward design objective is to produce a desirable response in tracking 
input commands when no random noises and disturbances are present. In particular, a 
fast and smooth tracking response during the transient phase when difficult maneuvers are 
being performed is the goal of the feedforward law. 

The main feedback design objectives are to suppress sensor noise, accommodate plant 
disturbances and provide stability in the presence of plant modeling uncertainties. 

Since the feedforward response must be fast while the feedback response must be 
relatively slow (to attenuate noise), combining these objectives into a single cost function 
produces conflicting demands. In this case, neither of the objectives is fully achieved. 
By separating the design objectives and decoupling the feedforward and feedback design 
processes, we are able to achieve both objectives fully. 

An important by-product of the SOFFT approach is that by removing the conflicting 
demands on the control system, we can design both feedforward and feedback control laws 
with less effort and in less time. The difficult and long process of fixing one problem only 
to find that another problem has popped up somewhere else is characteristic of conflicting 
demands. This effort seems to be mostly eliminated in the SOFFT approach. 
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Our main conclusion is that the concept of separating the feedforward and feedback 
objectives and decoupling the two designs as embodied in SOFFT works well. It pro- 
duces both fast and smooth command tracking as well as noise attenuation, stability and 
robustness. The design process also takes less effort. 

The SOFFT methodology can be used to design constant-gain feedforward and feed- 
back systems as well as variable-gain control systems. Our experience with constant-gain 
design shows that we can achieve perfect tracking when the plant matches the model used 
in the feedforward design. When the plant moves to a non-matching operating point, the 
response seems highly robust. 

A new measure of command tracking performance was developed. By analyzing these 
Z-plots at off-nominal operating points, we can predict the sensitivity or robustness of 
the total system in tracking commands. Z-plots provide an important tool for designing 
robust control systems. 

The Variable-Gain SOFFT methodology was used to design a flight control system for 
the F / A-18 aircraft. The purpose of using a Variable-Gain SOFFT controller was to extend 
the operating regime of the aircraft and to provide greater performance (flying/handling 
qualities) throughout the extended regime. The angle-of-attack (cc) command system 
designed brings the aircraft to within 3 degrees of its physical (actuator) limit. A detailed 
nonlinear simulation of the aircraft and control system displays excellent command tracking 
performance. 

Finally, all the gain computation algorithms, and the analysis tools shown in this 
report were obtained using our Computer-Aided Control System Design (CACSD) software 
tool ACET™. A new module to ACET has been developed in conjunction with the 
theoretical methodology development. ICS plans to market this new SOFFT design module 
in the coming months. 
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In conclusion, the high degree of success of the SOFFT design methodology has been 
rewarding. We recommend testing a Variable-Gain SOFFT flight control system in flight 
as a demonstration of this new technology. 
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APPENDIX A 


Multi-Model Robust SOFFT 

In section II, the SOFFT approach considered the control design problem for a single 
linear time-invariant (LTI) plant. Since most physical plants are nonlinear, it is more 
realistic to deal with a collection of LTI plants corresponding to a number of operating 
points. This collection of LTI plants may be referred as a multi-model plant and may 
also include plant modeling uncertainties such as uncertainties in parameter values or 
unmodeled dynamics. 

The objective in investigating multi-model plants is to design robust SOFFT con- 
trollers which have a satisfactory performance with all of the plant models. Thus, the 
controller deals with model uncertainties and operating point variations without changing 
any controller parameter values. It is important to note that the Variable-Gain SOFFT 
Controller also accommodates variations in operating points; however, the controller pa- 
rameters themselves vary to adapt to the plant variability. 

Consider a multi-model plant with M LTI models described by the system equations 

Xik+1 = X\ k "b F % Uik i 1 — * — M , (A.l) 


Uik — Cj Xik 


1 < t < M 


(A.2) 


In the single model case (i.e., when M = l), the SOFFT feedforward controller is 
obtained by solving an optimal control problem. Note that the SOFFT feedforward control 
system shown in Figure A.l is a dynamic compensator, not simply a gain matrix. 

Now, let the command model be given by 
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Vxk — C x x k (A. 5) 

with an input vector, u* A , of the form 

u* zk = -Klx* k - K z z k - K u u zk (A. 6) 

so that when u* xk is the input of the i tA plant model in (A.l) - (A.2), the resulting tracking 
error given by 


£\k — H yi t/ilc H z i Z k — H It Xik S z { z k (A. 7) 

is small for all the models, 1 < i < M. 

In this general formulation, the optimization is not only over the gain matrices K*,K Z 
and K u . We also want to find the dynamic compensator in (A.4) ; i.e., <f>* x and T* are also 
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psjt of the optimization although this was not originally proposed. This is a desirable 
feature, because a multi-model plant does not always identify (or favor) one particular 
model over the others. 

A second point of importance is that this formulation does not restrict the order (say 
Jij*) of the feedforward dynamic compensator. For example, the multi-model plant order 
may be 15, whereas we would like to design a feedforward controller of lower order, say 
4. This feature was also not originally proposed, but can be very important in practical 
applications, as has been stressed by Aaron Ostroff, of NASA, Langley Research Center. 
Finally, the lower order capability also provides a model reduction technique which should 
be investigated in the future. 

The cost function, J, may be selected as 


where 


N 


M 


J = 


2 (N + 1) £ 


£ fiE\\e ik \\l + rEWelWl 


k = 0 U=1 


(A.8) 


e k Hy y xk H z z k — H x x k U z z k (A. 9) 

The scalar weights /, prioritize the multiple models. Note that all the models are 
included in the cost function so that the design will try to reduce the tracking error for 
each model. More terms can be added to the cost function to achieve various additional 
goals. In particular, for order reduction, a slightly different cost is used. Here, we assume 
that the cost is quadratic. 

To develop an algorithm which determines the optimal controller parameters, we will 
formulate the general problem. We approach the problem in two steps. First, we consider 
the case in which u zk — 0. Then, we introduce a non-zero u zk . Now, consider the following 
open-loop feedforward model. 

2-Jfc+l = $0 % k "b u z k “b r oc u ok (A. 10) 
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where x* k is the feedforward controller state with n z - components, u* xk is the feedforward 
control vector with n ux components and u* ok is a feedforward compensator optimization 
control vector with n* 0 components. It should be noted that the purpose of u* ok is to 
optimize the feedforward dynamics; i.e., <j>* x and T*. The open-loop feedforward controller 
matrices <j> 0 and T ox may be selected as the nominal plant parameters when possible. In 
reduced-order control laws, they may be selected heuristically. 

The control vectors are of the form 


u' xk = -Klx\-K z z k 


(A.11) 


U ok — K ox Xu K oz Z)c 


(A.12) 


where the white noise u zk has been neglected. 

Suppose that the feedforward control, u* xk , is input to each of the plant models in a 
parallel fashion. We obtain the following augmented system. 
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From (A. 11) and (A.12), we select the output or feedback vector 


(A.13) 


y ‘ = fiM£) = (° o ::: l o ^ 
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(A.15) 


K* 


Kl K z \ 

k; x k oz ) 


(A.16) 


For the case of the feedforward design, it is reasonable to set the measurement noise 
covariance to zero. However, we have included these terms for generality. 

The augmented system in (A. 13) and (A. 14), coupled with the control constraint 
(A.15) and the cost function (A. 8) form an optimal output feedback control problem treated 
by the Principal Investigator previously. Therefore, this algorithm may now be used to 
obtain the augmented gain matrix, K*, in (A.16). 

Note that the cost function in (A. 8) is of quadratic form and may be computed as 
follows. 


H xi = H yi Ci 


1 < t < M 


(A.17) 


TT * TT* S~1* 

-*■* x ^ x 


(A.18) 


[H xl 0 

0 h x2 


H = 


V 0 


ffJ 0 

0 fj 


F = 


V 0 
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HxM 0 
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—Hzi \ 
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—H 2 m 
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0 J 

o\ 


f*I 


0 / 


(A.19) 


(A. 20) 


Q* = H T FH + Ql 


(A.21) 
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J = 5(jvTT) \ jt X ? «' x i + < T R ' < } (A.22) 

Consider the case in which u* k is not present. When the optimal K* is computed, the 
gains K * and K z determine u* xk in the feedforward control law in (A.4) and (A. 5) with 

<f>*x = <t>o , T* x = T ox . (A. 23) 

In the more general case, solving the optimization problem produces the gain K * and 
its partition in (A. 16). In this case, some manipulation yields 

<t>l = 4>o- r oc {k; x - lz K* x ) (a. 24) 

r* = r ox + r oc7z (a. 25) 

lz = K oz Kj [K z K t z \ ~ x (A. 26) 

For (A. 24) - (A. 26) to be valid, it is sufficient that 

r„ c {/ - K T Z [K z Jfj] _1 K z ) = 0 ■ (A. 27) 

While (A.27) will not hold in general, it will be satisfied when K z is non-singular. 

To include the effect of u z k , we augment the system given in (A. 13) by 

u zk+l = Ou zk + u UMk (A. 28) 

and couple it to z k through the command model (A. 3). Using previously developed results* 
to this augmented system, we obtain 

* A Combined Stochastic Feedforward and Feedback Control Design Methodology With Ap- 
plication to Autoland Design , NASA CR-4078, July 1987, p. 14. 
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k u = P x 1 r 5 Pxz r, w u {w u + vj- 1 

(A. 29) 


Px=T T x Px Tx + R* 

(A. 30) 


W n = E(u zk u T zk ) 

(A. 31) 

l 

i 

i 

E 

p _ ( Px Pxz \ 
\P,X Pzz) 

(A. 32) 

\ 


; I 

where P is optimal the cost matrix for the system in (A.13) and V u is the covariance of the | 

1 | 

measurement noise for u zk • Note that when V u = 0, which is usually valid in feedforward 

1 I 

control problems, W u cancels from (A. 29). ; 
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